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Oh! Abstract 
D , 

• We demonstrate that, together with the available experimental infor- 
mation, chiral symmetry determines the low energy behaviour of the tttt 

^ \ scattering amplitude to within very small uncertainties. In particular, the 

• threshold parameters of the S-, P-, D- and -F-waves are predicted, as 



well as the mass and width of the p and of the broad bump in the S'-wave. 
The implications for the coupling constants that occur in the effective La- 
grangian beyond leading order and also show up in other processes, are 
discussed. Also, we analyze the dependence of various observables on the 
mass of the two lightest quarks in some detail, in view of the extrapo- 
lations required to reach the small physical masses on the lattice. The 
analysis relies on the standard hypothesis, according to which the quark 
condensate is the leading order parameter of the spontaneously broken 
symmetry. Our results provide the basis for an experimental test of this 
hypothesis, in particular in the framework of the ongoing DIRAC experi- 
ment: The prediction for the lifetime of the ground state of a tt+tt" atom 
reads r = (2.90 ± 0.10) 10"^^ sec. 

Pacs: ll.30.Rd, 12.38.Aw, 12.39.Fe, 13.75.Lb 

Keywords: Roy equations, Meson-meson interactions, Pion-pion scattering, 
Chiral symmetries 
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1 Introduction 

The study of tttt scattering is a classical subject in the field of strong interactions. 
The properties of the pions are intimately related to an approximate symmetry 
of QCD. In the chiral limit, where and vanish, this symmetry becomes 
exact, the Lagrangian being invariant under the group SU(2)axSU(2)L of chiral 
rotations. The symmetry is spontaneously broken to the isospin subgroup SU(2)v. 
The pions represent the corresponding Goldstone bosons. 

In reality, the quarks are not massless. The theory only possesses an ap- 
proximate chiral symmetry, because ruu and rrid happen to be very small. The 
consequences of the fact that the symmetry breaking is small may be worked out 
by means of an effective field theory The various quantities of interest are 
expanded in powers of the momenta and quark masses. In the case of the pion 
mass, for instance, the expansion starts with 

M^+ = {mu + md)B + Oim^), 5 = |(0|mm |0)| , (1.1) 

where F is the value of the pion decay constant in the chiral limit, mu,mii 0. 
The formula shows that the square of the pion mass is proportional to the product 
of rUu + md with the order parameter (OImm |0). The two factors represent quan- 
titative measures for explicit and spontaneous symmetry breaking, respectively. 
If the explicit symmetry breaking is turned off, the pions do become massless, as 
they should: The symmetry is then exact, so that the spectrum contains three 
massless Goldstone bosons, while all other levels form massive, degenerate isospin 
multiplets. 

The properties of the Goldstone bosons are strongly constrained by chiral 
symmetry: In the chiral limit, the scattering amplitude vanishes when the mo- 
menta of the pions tend to zero. To first order in the symmetry breaking, the 
S'-wave scattering lengths are proportional to the square of the pion mass 0: 

where a| stands for the scattering length in the isospin / channel with angular 
momentum £. The two low energy theorems ( p..2|) are valid only at leading order 
in a series expansion in powers of the quark masses: The next-to-leading order 
corrections were calculated in 0, and even the next-to-next-to-leading order 
corrections are now known lEI. 
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In the following, we exploit the fact that analyticity, unitarity and crossing 
symmetry impose further constraints on the scattering amplitude. These were 
analyzed in detail in on the basis of the Roy equations and of the exper- 
imental data available at intermediate energies. The upshot of that analysis is 
that Og and are the essential low energy parameters: Once these are known, the 
available experimental data determine the low energy behaviour of the vrvr scat- 
tering amplitude to within remarkably small uncertainties. As discussed above, 
chiral symmetry predicts exactly these two parameters. Hence the low energy be- 
haviour of the scattering amplitude is fully determined by the experimental data 
in the intermediate energy region and the theoretical properties just mentioned: 
analyticity, unitarity, crossing symmetry and chiral symmetry. 

The resulting predictions for the S-wave scattering lengths were presented 
already 0. The purpose of the present paper is to (i) discuss the analysis that 
underlies these predictions in more detail, (ii) present the results for the threshold 
parameters of the P , and F-waves, (iii) give an explicit representation 
for the S- and P-wave phase shifts and (iv) extract the information about the 
coupling constants of the effective Lagrangian. 

Several authors have performed a comparison of the chiral perturbation the- 
ory predictions with the data, in particular also in view of a determination of the 
effective coupling constants ii and £2 lH^liU- Stern and collaborators ^ 



investigate the problem from a different point of view, referred to as "Generalized 
Chiral Perturbation Theory" . These authors treat the S'-wave scattering lengths 
as free parameters and investigate the possibility that their values strongly de- 
viate from those predicted by Weinberg. In the language of the effective chiral 
Lagrangian, this scenario would arise if the standard estimates for the effective 
coupling constant £3 were entirely wrong: The quark condensate would then fail 
to represent the leading order parameter of the spontaneously broken chiral sym- 
metry. Indeed, these estimates rely on a theoretical picture that has not been 
tested experimentally. 

On the experimental side, the situation is the following. As shown in early 
numerical analyses of the Roy equations ||25|, only data sufficiently close to thresh- 
old can provide significant bounds on the scattering lengths. The often quoted 
values a° = 0.26 ± 0.05, = -0.028 ± 0.012 p|, ^ mainly rely on the 3 • 10^ 
K vrvrez/ decays collected by the Geneva-Saclay collaboration, which provided 



its final results in 1977 |2^. There are new data from Brookhaven [29, 30], where 
more than 4 • 10^ K^^ decays are being analyzed, and the low energy behaviour of 



the relevant form factors is now also known much better p^ , |22| . As will be dis- 
cussed in section |T3|, the preliminary results of this experiment indeed reduce the 
uncertainties significantly. A similar experiment is proposed by the NA48 collab- 
oration at CERN [Q. Unfortunately, the data taking at the DA$NE facility is 
delayed, due to technical problems with the accelerator. A beautiful experiment 
is under way at CERN |3^, which is based on the fact that tt+tt^ atoms decay 
into a pair of neutral pions, through the strong transition tt+tt^ tt^tt^. Since 
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the momentum transfer nearly vanishes, only the scattering lengths are relevant: 
At leading order in isospin breaking, the transition amplitude is proportional to 
Cq — Oq. The corrections at next-to-leading order are now also known |^, as 
a result of which a measurement of the lifetime of a ir'^n" atom amounts to a 
measurement of this combination of scattering lengths. Finally, we mention the 
new data on pion production off nucleons, obtained by the CHAOS collabora- 
tion at Triumf ||36|. The scattering lengths may be extracted from these data 
by means of a Chew-Low extrapolation procedure. Chiral symmetry, however, 
suppresses the one-pion exchange contribution with a factor of t, so that a careful 
data selection is required to arrive at a coherent Chew-Low fit. It yet remains to 
be seen whether these data permit a significant reduction of the uncertainties in 
the experimental determination of aj] and a^. 

The experiments mentioned above are of particular interest, because they offer 
a test of the hypothesis that the quark condensate represents the leading order 
parameter of the spontaneously broken symmetry: If the predictions obtained 
in the present paper should turn out to be in contradiction with the outcome 
of these experiments, the commonly accepted theoretical picture would require 
thorough revision. 



2 Chiral representation 

Throughout the present paper we work in the isospin limit: We disregard the 
e.m. interaction and set m„ = = m. The various elastic reactions among two 
pions may then be represented by a single scattering amplitude A{s,t,u). Only 
two of the Mandelstam variables are independent, s + t + u = 4M^ and, as a 
consequence of Bose statistics, the amplitude is invariant under an interchange 
of t and u. 

As discussed in detail in ref . , chiral perturbation theory allows one to study 
the properties of the vrvr scattering amplitude that follow from the occurrence 
of a spontaneously broken approximate symmetry. The method is based on a 
systematic expansion in powers of the momenta and of the light quark masses. 
We refer to this as the chiral expansion and use the standard bookkeeping, which 
counts the quark masses like two powers of momentum, m = 0{p'^). 

The two loop representation of the scattering amplitude given in yields the 
first three terms in the chiral expansion of the partial waves: 

t[(s) = i{s)2 + t[{s), + t{{s), + 0(/) (2.1) 

At leading order, only the S- and P-waves are different from zero: 
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In the low energy expansion, inelastic reactions start showing up only at 0{p^). 
The unitarity condition therefore reads: 



Imtiis) = a{s) \tl{s)\' + 0{/) , a(s) = ^1 - — ^ . (2.3) 

The condition immediately implies that the imaginary parts of the two loop 
amplitude may be worked out from the one-loop representation: 

lmtl{s)^a{s)tl{s)2 {tl{s)2 + 2Retl{s),} + Oip') (2.4) 

The formula shows that, at low energies, the imaginary parts of the partial waves 
with i > 2 are of order and hence beyond the accuracy of the two loop 
calculation. 

Stated differently, the imaginary part of the two loop representation is due 
exclusively to the S- and P-waves. This implies that, up to and including O(p^), 
the chiral representation of the scattering amplitude only involves three functions 
of a single variable: 

A{s, t, u) = C{s, t, u) + 327r {i U\s) + | (s - m) U\t) + | (s - t) U\u) 

+ 1 u\t) + \ U\u) - \ U\s)} + 0(/) . (2.5) 

The first term is a crossing symmetric polynomial in s, rt, 

C (s, t, ti) = Ci + S C2 + C3 + {t — uY C4^ + + S {t — uY Cq . (2.6) 

The functions U^{s), U^{s) and C/^(s) describe the "unitarity corrections" associ- 
ated with s-channel isospin 7 = 0, 1, 2, respectively. In view of the fact that the 
chiral perturbation theory representation for the imaginary parts of the partial 
waves grows with the power \mt{{s)Q oc s'^, we need to apply several subtractions 
for the dispersive representation of these functions to converge. It is convenient 
to subtract at s = and to write the dispersion integrals in the form 



U\s) = 




POO 




/ ds' 


TT , 


JAMl 






roo 


U\s) = 




/ ds' 


TT . 


JAMl 






roo 


U\s) = 




/ ds' 


TT . 





s'\s' - 4M2) {s' - s) 
, <7(.')^^(«')2 0g(O2 + 2Retg(A} 



(2.7) 



s>\s' - s) 



The subtraction constants are collected in the polynomial C(s, t,M). Alterna- 
tively, we could set C(s, li) = and book the subtraction terms as polynomial 
contributions to V^i^s). U^{s). U'^{s). The decomposition of C{s,t,u) into a set 
of three polynomials of a single variable is not unique, however, so that we would 
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have to adopt a convention for this sphtting - we find it more convenient to work 
with the above representation of the amphtude. 

The specific structure of the unitarity correction given above was noted al- 
ready in . It is straightforward to check that the exphcit resuh of the full two 
loop calculation described in ^ is indeed of this structure. The essential result 
of that calculation is the expression for the polynomial part of the amplitude, 
in terms of the effective coupling constants. The corresponding formulae, which 
specify how the coefficients Ci, . . . , Cg depend on the quark masses, are given in 
appendix 0. These, in particular contain Weinberg's low energy theorem, which 
in this language states that the expansion of the coefficients ci and C2 starts with 



The two loop calculation specifies the expansion of these two coefficients up to 
and including next-to-next-to-leading order. 

3 Phenomenological representation 

As shown by Roy [0, the fixed-t dispersion relations for the isospin amplitudes 
can be written in such a form that they express the vrvr scattering amplitude 
in terms of the imaginary parts in the physical region of the s-channel. The 
resulting representation for A{s,t,u) contains two subtraction constants, which 
may be identified with the scattering lengths and Oq. Unitarity converts this 
representation into a set of coupled integral equations, which we recently exam- 
ined in great detail 0. In the present context, the main result of interest is that 
the representation allows us to determine the imaginary parts of the scattering 
amplitude in terms of Oq and Oq. Since the resulting representation is based on 
the available experimental information, we refer to it as the phenomenological 
representation. 

In the following, we treat the imaginary parts of the partial wave amplitudes 
as if they were completely known from phenomenology - we will discuss the 
uncertainties in these quantities as well as their dependence on Og and Oq in 
detail, once we have identified the manner in which they enter our predictions 
for the scattering lengths. 

The chiral representation shows that the singularities generated by the imag- 
inary parts of the partial waves with i > 2 start manifesting themselves only 
at 0{p^). Accordingly, we may expand the corresponding contributions to the 
dispersion integrals into a Taylor series of the momenta. The singularities due to 
the imaginary parts of the S- and P-waves, on the other hand, start manifesting 
themselves already at O(p^) - these cannot be replaced by a polynomial. The 
corresponding contributions to the amplitude are of the same structure as the 
unitarity corrections and also involve three functions of a single variable. It is 
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convenient to subtract the relevant dispersion integrals in the same manner as 
for the chiral representation: 





TT 


/ ds' 

JAMl 


W\s) = 


53 
TT 


rco 

/ ds' 


W\s) = 


TT 


fOO 

/ ds' 



, lmtl{s') 

(3.1) 



s'\s' - s) ' 

IlllfJ(,' 



Imtlis') 



Since all other contributions can be replaced by a polynomial, the phenomeno- 
logical amplitude takes the form 

A{s, t, u) = IQiral + ^ (2a° - Sa^) s + P{s, t, u) (3.2) 

+327r {|W^°(s) + |(s - ^/)fy^(i) + |(s - t)fy^(^/) 
+ |W^2(t) + |f^2(«)-iW^2(.)} + 0(/). 

We have explicitly displayed the contributions from the subtraction constants Qq 
and Cq. The term P{s, t, u) is a crossing symmetry polynomial 

P(s, t, u) = pi +p2 S +p3 +p4 (t - uf +P5 S^ -\-pQ s{t - ■ (3.3) 

As demonstrated in the appendix, its coefficients can be expressed in terms of 
the following integrals over the imaginary parts of the partial waves: 

ji _ ^ (2^ + 1) /~ , Imti(g) 

1 /-^ , 2 Im tO(s) +4 Im t2(s) 
' (is 



= V (2/ + !)£(£+!)- / 



^=2 

The explicit expressions read 



4M2 9s3(s-4M2) 



Pi = -1287rM^^ (7J + 72 + 2Mll\ + 2M2/2 + SM,^/^) , 



p3^Y (^^0 ~ ~ ~ leM^J? - 42M2/J + 22M2/2 - 72M^il + 2AMIh) , 
P4 = Stt (/J + 7^ + 2Mll\ + 2M^72 _ 24717^7^^ , (3.5) 

p^ = ^ (87° + 97} - 117? - 32M27° + AAMlll - 677) , 



3 

P6 = 47r (/} - 3/2 + l2Mlll + 277) . 



The fact that, at low energies, the scattering amphtude may be represented in 
terms of integrals over the imaginary parts that can be evaluated phenomeno- 
logically, was noted earlier, by Stern and collaborators |2^. These authors also 
worked out the implications for the threshold parameters and the effective cou- 
pling constants of the chiral Lagrangian and we will compare their results with 
ours, but we first need to specify the framework we are using. 



4 Matching conditions 

In the preceding sections, we have set up two different representations of the 
scattering amplitude: One based on chiral perturbation theory and one relying 
on the Roy equations. The purpose of the present section is to show that, in 
their common domain of validity, the two representations agree, provided the 
parameters occurring therein are properly matched. 

The chiral and phenomenological representations are of the same structure. 
The coefficients of the polynomials C{s, t, u) and P(s, t, u) are defined differently 
and, instead of the functions U^{s) occurring in the chiral representation, the 
phenomenological one involves the functions W\s). The latter are defined in 
eq. (|3.1| ), as integrals over the imaginary parts of the physical 5*- and P-waves. 

The key observation is that, in the integrals ( |3.1| ), only the region where s' is 
of order matters for the comparison of the two representations. The remainder 
generates contributions to the amplitude that are most of order p^. Moreover, for 
small values of s', the quantities Imt^(s') are given by the chiral representation in 
eq. ( |2.4D except for contributions that again only manifest themselves at 0{p^). 
This implies that the differences between the functions W\s) and U\s) are 
beyond the accuracy of the chiral representation: 

= f/°(s) + 0{p^) , 
W\s) = U\s) + Oip^) , (4.1) 
W\s) = U\s) + 0{p^). 

Hence the two representations agree if and only if the polynomial parts do, 

4:71 — 

C{s, t, u) = IQiial + ^ (2a0 - 5a^) s + P{s, t, u) + 0{p') . 
This implies that the coefficients of C{s,t,u) and P{s,t,u) are related by 

4:71 

Cl = 167Tal +pi + 0(/) , C2 = ^ (2a[] - Sa^) +p2 + 0(/) , 

C3 = P3 + 0{p^), c, = p, + 0{p^), (4-2) 

C5 = P5 + , Cq = Pq + 0{p'^) . 
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The chiral representation specifies the coefficients ci, . . . , cq in terms of the 
effective couphng constants, while the quantities pi, are experimentally 
accessible. Since the main uncertainties in the latter arise from the poorly known 
values of the scattering lengths Oq, Oq, the above relations essentially determine 
the coefficients ci, . . . , cg in terms of these two parameters. 

5 Symmetry breaking in the effective Lagrangian 

As discussed in section 0, unitarity fully determines the scattering amplitude to 
third order of the chiral expansion, in terms of the couplings constants occurring 
in the derivative expansion of the effective Lagrangian to 0{p^), 

Ceff = C2 + U + + ... (5.1) 

The leading term C2 only contains F and = 2mB. The vertices relevant for 
HTC scattering involve the coupling constants £1, £2, ^3, ^4 from £4, and Cq generates 

6 further couplings: ri, . . . , rg. We need to distinguish two different categories of 
coupling constants: 

a. Terms that survive in the chiral limit. Four of the coupling constants that 

enter the two loop representation of the scattering amplitude belong to this 
category: ii, £2, r^, Tq. 

b. Symmetry breaking terms. The corresponding vertices are proportional to a 

power of the quark mass and involve the coupling constants £3, £4, ri, r2. 

The constants of the first category show up in the momentum dependence of the 
scattering amplitude, so that these couplings may be determined phenomenolog- 
ically. The symmetry breaking terms, on the other hand, specify the dependence 
of the amplitude on the quark masses. Since these cannot be varied experimen- 
tally, information concerning the second category of coupling constants can only 
be obtained from sources other than vrvr scattering. In part, we are relying on 
theoretical estimates here. Although these are rather crude, the uncertainties do 
not significantly affect our results, for the following reason. 

The quark masses rUu, m^, which are responsible for the symmetry breaking 
effects, are very small compared to the intrinsic scale A of the theory, which is of 
order 500 MeV or 1 GeV. The group SU(2)r xSU(2)l therefore represents a nearly 
perfect symmetry of the QCD Hamiltonian. In the isospin limit, the symmetry 
breaking effects are controlled by the ratio m/A, with m = ^ (m„ + nid). In view 
of m ~ 5 MeV, the expansion parameter is of the order of 10~^, indicating that 
the expansion converges very rapidly. 

In the framework of the effective theory, it is convenient to replace powers 
of m by powers of and to identify the intrinsic scale A with AirF.^. The 
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expansion parameter m/A is then replaced by 



The numerical valueQ ^ = 1.445 ■ 10~^ confirms the estimate just given. 

We know of only one mechanism that can upset the above crude order of 
magnitude estimate for the symmetry breaking effects: The perturbations gener- 
ated by the quark mass term in the QCD Hamiltonian, m^fZw + m^Jc?, may be 
enhanced by small energy denominators. Indeed, small energy denominators do 
occur: 

(i) In the chiral limit, the pions are massless, so that the straightforward 
expansion in powers of the quark masses leads to infrared singularities. For a 
finite pion mass, these singularities are cut off at a scale of the order of M^- and 
the divergences are converted to finite expressions that involve the logarithm of 
M^-. The most important contributions of this type are generated by the vertices 
contained in the leading order effective Lagrangian, which are fully determined 
by and M^.. Accordingly, the coefficients of the leading chiral logarithms do 
not involve any unknown constants. In those cases where this coefficient happens 
to be large, the symmetry breaking effects are indeed enhanced, so that the above 
rule of thumb estimate then fails. 

(ii) States that remain massive in the chiral limit may give rise to small energy 
denominators if their mass happens to be small. In the framework of chiral 
perturbation theory, the occurrence of such states manifests itself only indirectly, 
through the fact that some of the effective coupling constants are comparatively 
large. The p-meson represents the most prominent example and it is well-known 
that some of the coupling constants (for instance ii and £2) are dominated by the 
contribution from this state In fact, for all of those effective couplings that 
have been determined experimentally, the observed magnitude is well accounted 
for by the hypothesis that they are dominated by the resonances seen at low 
energies 



6 Low energy theorems 

As the two loop formulae are rather lengthy, we first discuss the principle used 
to arrive at the prediction for the S'-wave scattering lengths at one loop level, 
where the algebra is quite simple. The first order corrections to the two low 
energy theorems ( p.8|) are readily obtained from the formulae given in appendix 
|B[ Expressed in terms of the scale invariant effective coupling constants ii, . . . ,£4, 

^ Throughout this paper, we identify with the mass of the charged pion and use F-^ = 
92.4 MeV @. 
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introduced in [Q, the result reads: 

If. . / 4 ^ ^ ^ 67 \ 



c. = ^(1+^(^-3^1 + 24-^; +0(^1 . (6.1) 

The corrections involve both types of couplings: £i is of type a. and can thus be 
determined from the momentum dependence of the scattering amplitude, while 
^3 and £4 are of type b.. Indeed, both ii and £2 show up in the terms proportional 
to and (t — m)^: 

1 ill I2 471 1 [£2 1271 

These formulae show that, up to and including terms of order ^, the quantities 
C,^Fl[c2 + AMl{c,-c,)}, C2^^{-ci+4M^^(c3-C4)} (6.2) 
exclusively contain the symmetry breaking couplings £3 and £4: 

C. = l + i[2h-^^] + 0{e), C72 = l+e|| + 2£4-y}+0(a. 

In the following, we analyze the low energy theorems for the S-wave scattering 
lengths by means of the quantities Ci and C2 defined in eq. ( |6.2| ). The one for 
2ag — Sag, for instance, is obtained by inserting the matching relations ( ^.2| ) in 
the definition of Ci and solving for the scattering lengths. The result reads 

2a° - Sag = f-^ + M'^ a, + 0(M«) , (6.3) 

where ai collects the contributions from the phenomenological moments, 

ai = IQ M^{8P^ + 9 II- 11 II - 36M^ II - Q H) . (6.4) 
The analogous low energy theorems for ag and read 

7M2 

ao = Co + ag + 0(M«) , (6.5) 

ag = -Y^C^2 + M,^a2 + 0(M«), 
where Co is a combination of Ci and C2, 

Co = ^(12Ci-5C2), (6.6) 
12 



while ao, a2 again stand for a collection of moments 

tto = ^ (5/o + 10/^ + 28M^ll + 2AMll\ - IGM^ll - 96M^ll + 6M^H) , 
a2 = ^ (/S + 2 72 - 4M270 - 6M^ll + lOM^ll + 24M^ll + 6M^H) . (6.7) 



The relations (|6.2| )- (|67^ ) specify the S-wave scattering lengths in terms of Ci, C2 
and the phenomenological moments 7^ and H . Note that these contain infrared 
singularities. Their chiral expansion starts with the contributions generated by 
the square of the tree level amplitudes: 

227 

+ 0(M-^), 



/? = 


101 

M^K 


+ 0(1), 


n 


I\ = 


2 

MIK 


+ 0(1), 


n 


n = 


14 


+ 0(1), 


2 


H = 


0(1), 




K 



lAMtK 



^ + 0{Mf) , (6.8) 



7M^K 
61440 vr^F^ . 

The evaluation of the moments requires phenomenological information. Since 
the behaviour of the imaginary parts near threshold is sensitive to the scattering 
lengths we are looking for, the same applies for these moments. In the narrow 
range of interest, the dependence is well described by the quadratic formulae in 
appendix 0, which yield 

M^ao = .0448 + .30 Aa° - .37 Aa^ + .5 (Aa[!)^ - 1.2 Aaj^Aa^ + 1.8 {Aalf 
M^ai = .0619 + .48 Aa° - .26 Aa^ + .86 (Aa")^ - 1.7 Aa[jAa^ + .3 {Aalf 
M^a2 = .00553 + .023 Aa° - .095 Aal - .1 Aa^Aa^ + .7 {AalY , (6.9) 

with AaO = a[] - 0.225, Aag = + 0.03706. 



7 The coupling constants £3 and £4 

The representation of the S-wave scattering lengths derived in the preceding 
section splits the correction to Weinberg's leading order formulae into two parts: 
a correction factor C„, which at first nonleading order only involves the coupling 
constants and ^4 and a term a„ that can be determined on phenomenological 
grounds. 

The significance of the coupling constants £3 and £4 is best seen in the ex- 
pansion of and in powers of the quark mass. The relation of Gell-Mann, 
Oakes and Renner states that the expansion of starts with a term linear 
in m. The coupling constant £3 determines the first order correction: 

= m2{1- ie4 + 0(e')}, M^ = 2Bm. (7.1) 
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The constant B stands for the value of | (0|fZM |0) in the chiral hmit. Note that 
£3 contains a chiral logarithm, £3 = — InM^ + 0(1). The coupling constant I4, 
which also contains a chiral logarithm with unit coefficient, £4 = — lnM^ + 0(l), 
is the analogous term in the expansion of the pion decay constant, 

F^ = F{l + Ch + 0{e)}, (7.2) 

where F is the value of F^^ in the chiral limit. 

The same two coupling constants also show up in the scalar form factor 

{n{p') I rriuuu + mddd\ tt{p)) = {1 + | {r^), t + 0{t^)} . (7.3) 

The value of the matrix element at t = is the pion a-term. According to the 
Feynman-Hellman theorem, it is given by a-jr = mdM^/dm. The relation ( |7.1| ) 
thus shows that £3 also determines the a-term to first nonleading order: 

a. = A4 {1 - (£3 - 1) + a. + 0{e)} . (7.4) 

Moreover, chiral symmetry implies that the same coupling constant that deter- 
mines the difference between F^^ and F also fixes the scalar radius at leading 
order of the chiral expansion |]^: 

(-^). = ^{^"4-§ + eA. + 0(a}. (7.5) 

We may therefore eliminate £4 in favour of the scalar radius and rewrite the 
correction factors in the form 



3 ' " 14 I 525 
C^i = l + ^(r\ + || + e^A, + 0(a, (7.6) 

C^2 = l + ^(r\ + |{£3-g} + e^A, + 0(a, 

with Aq = (12 Ai — 5 A2)/7. The first order corrections are then determined by 
{r'^)s and £3, while Aq, Ai and A2 represent the two loop contributions. The 
scalar form factor is also known to two loops The explicit expressions for 
the second order corrections are given in appendix 0. 

For the numerical value of the scalar radius, we rely on the dispersive eval- 
uation of the scalar form factor described in ref. We have repeated that 
calculation with the information about the phase shift Sq{s) obtained in ref. ||^. 
In view of the strong final state interaction in the S'-wave, the scalar radius is 
significantly larger than the electromagnetic one, (r 

2)_ = 0.439 ± 0.008 fm^ g. 

The result reads 

(r^), = 0.61 ±0.04 fm% (7.7) 
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where the error is our estimate of the uncertainties to be attached to the dispersive 
calculation. The number confirms the value given in ref. 0] and is consistent 



with earlier estimates of the low energy constant £4, based on the symmetry 



breaking seen in Fk/F^^ or on the decay K — >■ tt^z/ |^2[, but is more accurate. 
It corresponds to {'f'^) s = 0.102 ± 0.007, so that the contribution from the 
scalar radius represents a correction of order 10%, in Co, Ci, as well as in C2. 

The crucial parameter that distinguishes the standard framework from the 
one proposed in ref. is £3. The value of this coupling constant is not known 
accurately. Numerically, however, a significant change in the prediction for the 
scattering lengths can only arise if the crude estimate 

£3 = 2.9 ± 2.4 (7.8) 

given in ref. should turn out to be entirely wrong: With this estimate, the 
contribution from £3 to and Oq is of order 0.002 and 0.001, respectively. We 
do not make an attempt at reducing the uncertainty in £3 within the standard 
framework, because it barely affects our final result. Instead, we will explicitly 
display the sensitivity of the outcome to this coupling constant. 



8 Results for and at one loop level 

We first drop the two loop corrections A„. Inserting the values (r'^)s = 0.61 fm^ 
and £3 = 2.9, the low energy theorems ( |7.6| ) yield 

Co = 1.092 , Ci = 1.103 , C2 = 1.117 . (8.1) 

The correction factor Ci is fully determined by the contribution from the scalar 
radius. The numerical values of Cq and C2 differ little from Ci: The estimate 
(|7.8|) implies that the contributions from the coupling constant £3 are very small, 
so that these terms are also dominated by the scalar radius. Inserting the values 
( |6.9| ), ( |8.1| ) in the relations ( |6.5| ) and solving for aQ,al, we then get 

a[! = 0.2195, a^ = -0.0446, 2ag - 5a^ = 0.662 . (8.2) 

These numbers are somewhat different from those obtained in [Q, which are 
also based on the one loop representation of the scattering amplitude. In fact, 
even if the two loop corrections A„ are dropped, the formulae ( |6.5D for the S- 
wave scattering lengths differ from those given in ref. 0]. In the case of aj], for 
example, the formula given there reads 

where and a"^ are the D-wave scattering lengths. As far as the contributions 
proportional to (r^)^ and £3 are concerned, the expression is the same, but instead 
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of the phenomenological moments contained in cto, the above formula contains 
the term 



tto < — ^ — (^2 + 2ai) + ^^^r^ . (8.3) 

^ V 2 ^ 2; ^ g^20 vr^F^ ^ ^ 

Indeed, the D-wave scattering lengths may be expressed in terms of moments, 
up to and including contributions of first nonleading order. Projecting the phe- 
nomenological representation (|3.2|) onto the D-waves, we find that the functions 
W\s) do not contribute to the scattering lengths, while the contribution from 
the background polynomial reads 

a° = ^ {n + 3/^ + 5/1 - AM'Ml - 3/i + 5/1) + 30M'^h} + 0{Mt) , (8.4) 
= A |2/1 - 3/1 + /I - AMl{2ll + 3/} + 7?) + 2AMIh} + 0{Mt) . 

The comparison with the exact representation for the /)-wave scattering lengths 
given in shows that the contributions from the imaginary parts of the S- and 
P-waves can be represented in terms of the moments and the coefficients agree 



with those above. The formula ( |8.4| ) includes the contributions from the higher 
partial waves, up to and including corrections of first nonleading order. In the 
difference, 

the leading moments cancel, but the terms with /{, and H remain: 
737 A/^ 

Aa° = - gy^Q^3^, +SMl (8/1 + 4/1 + 4/1 - //) - 128M«/1 . (8.5) 

The low energy expansion of the moments of eq. (|6.8|) shows that the contributions 
of O(M^) in Aoq indeed cancel out, demonstrating that the formula given in 
ref. [§| agrees with our representation, up to terms that are beyond the algebraic 
accuracy of that formula. 

Numerically, however, the leading order terms represent a rather poor approx- 
imation for the moments, so that there is a numerical difference: The numerical 
values of the moments are given in appendix ^. Inserting these in (|8.4|) , we obtain 
al = 1.76- 10"=^ M"^, al = 0.171-10-3 M-^ so that the one loop formula of ref . 
yields = 0.205, instead of the value = 0.2195 given above. The difference 
arises because we are matching the chiral and phenomenological representations 
differently: We represent the amplitude in terms of three functions of a single 
variable s and match the coefficients of the Taylor expansion at s = 0. In ref. 0], 
the one loop formulae for the various scattering lengths were obtained by directly 
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evaluating the chiral representation at threshold - in other words, the matching 
was performed at s = 4M^ rather than at s = 0. 

We emphasize that the above discussion in the framework of the one loop ap- 
proximation only serves to explicitly demonstrate that the choice of the matching 
conditions is not irrelevant. Admittedly, in our final analysis, where we will be 
working at two loop accuracy, the noise due to that choice is significantly smaller. 

9 Infrared singularities 

From a purely algebraic point of view, the manner in which the matching is 
done is irrelevant, as long as it is performed in the common region of validity of 
the chiral and phenomenological representations. We could also match the two 
loop representation to the phenomenological one at threshold and would then 
obtain a formula analogous to the one given in 0], but now valid to next-to- 
next-to-leading order. Alternatively, we could match the two representations of 
the scattering amplitude at the center of the Mandelstam triangle - the result 
would only differ by contributions that are beyond the accuracy of the chiral 
representation. 

There is a good reason for preferring the procedure specified above to a match- 
ing at threshold: The branch cut required by unitarity starts there. The mod- 
ifications of the tree level result generated by the higher order effects are quite 
large at threshold, because they are enhanced by a small energy denominator. 
Indeed, contains a chiral logarithm with an unusually large coefficient: 



The phenomenon gives rise to an exceptionally large correction that violates the 
rule of thumb of section |^ by an order of magnitude: The one-loop correction 
increases the tree level prediction by about 25% ! 

At the center of the Mandelstam triangle, the amplitude also contains a chiral 



The coefficient is less than half as big as the one in af], but it still represents a 
sizeable correction. 

In our matching procedure, we replace Oq and Og by Cq and C2 and at the 
same time also eliminate £4 in favour of the scalar radius. What matters for the 
convergence properties of the quantities appearing in our matching conditions are 
the infrared singularities contained in 




logarithm (sq = |M^): 





kf^ 1 
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The coefficients occurring here are remarkably small. The term Ci — ■|M^(r^)^ 
does not contain a chiral logarithm at all. We can therefore expect that, for 
the quantities that are relevant for the determination of the S'-wave scattering 
lengths, the perturbation series converges very rapidly, much more so than for a 
matching at threshold or at the center of the Mandelstam triangle. As we will 
see, this is indeed born out by the numerical analysis. 



10 Estimates for symmetry breaking at 0{p^) 

We now extend the analysis to next-to-next-to-leading order. For that purpose, 
we need an estimate for the symmetry breaking couplings ri ... , r4 and of 
Cq, which enter the low energy theorems for Cq,Ci,C2 at order M^, as well as 
the relation between the scalar radius and the coupling constant £4. The corre- 
sponding correction terms Ao,Ai,A2 are listed in ( |C.2| ). In the normalization 



used there, the resonance estimates of refs. 0, 21, 431 amount to 



fi~-1.5, ^2^3.2, f3~-4.2, f4~-2.5, fs2--0-7- (10.1) 
Inserting these numbers in (|C.2|) , we obtain a shift in Co, Ci, C2 by —0.3, —0.5 and 



—0.8 permille, respectively. This confirms the expectation that the effects due to 
the symmetry breaking coupling constants r„ are tiny. Since the scale is set by 
the scalar or pseudoscalar non-Goldstone states contributing to the relevant sum 
rules, Ms ~ 1 GeV, the corresponding corrections are of order M^/Mj ~ 4 ■ 10~^. 
In the SU(2) framework we are using here, the KK continuum also contributes 
to the effective coupling constants, but in view of AM"^ ~ M^, the corresponding 
scale is even somewhat larger. In the following, we assume that the estimates in 
equation ( |10.1 ) are valid to within a factor of two. 



In the case of ri, . . . , r4, the main uncertainty stems from the vrvr continuum 
underneath the resonances, that is from the chiral logarithms. Since the formulae 
(|C.2| ) are quadratic in these, the scale dependence of those coupling constants is 
rather pronounced. This can be seen by varying the scale fi, at which the running 
coupling constants are assumed to be saturated by the resonance contributions. 
For 0.5 GeV < fi < 1 GeV, the corrections vary in the range 

0.002 < e'Ao ^ 0.005, -0.001 < ^'Ai ^ 0.003, -0.005 < e'A2 < 0.001. 

In the representation { \l.5\j for the scalar radius, the two loop correction 
represents an effect of first order. Estimating the magnitude in the same man- 
ner as for Ao,Ai,A2, the result varies in the range 0.18 < ^A^ ^ 0.28. The 
correction thus shifts the scalar radius by 0.04 ± 0.01 fm^. 

In the following, the central values are calculated by using the resonance 
estimates ( |10.1| ) at the scale /i = Mp. For some of the quantities analyzed in 
the present paper, the result is insensitive to the uncertainties inherent in these 
estimates, but in some cases, they even dominate our error bars - we will discuss 
the sensitivity of the various results in detail. 
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0.220 


-0.0444 


-0.36 


4.31 


4.39 


3.8 


1.0 




0.002 


0.0003 


0.04 


0.02 


0.19 


0.05 


0.03 




0.004 


0.0009 


0.01 


0.00 


0.02 


0.01 


0.00 




0.001 


0.0002 


0.51 


0.10 


0.10 


1.04 


0.10 


exp 


0.001 


0.0002 


0.29 


0.04 


0.03 


0.12 


0.02 


tot 


0.005 


0.0010 


0.59 


0.11 


0.22 


1.05 


0.11 



Table 1: Solution of the matching conditions. The first row contains the central 
values. The next four rows indicate the uncertainties in this result, arising from 
the one in (r^)^, £3, r„ and in the experimental input used in the Roy equations. 
The last row is obtained by adding these up in quadrature. 

11 Final results for Oq and 

We are now in a position to describe the determination of Oq and Cq at two loop 
accuracy. Our matching conditions identify two different representations for the 



coefficients ci, ... ce : the chiral representation specified in equation ( |B.2|) and 
the phenomenological one in ( |4.2| ). For the evaluation of the S'-wave scattering 
lengths, only the first four coefficients are relevant. For these, the chiral represen- 
tation involves the effective coupling constants ii, £2, £3, ^i, ^i, r2, r^, r4, while 
the phenomenological representation contains only the two parameters Oq and Oq, 
which enter explicitly as well as implicitly, through the moments pi, . . . ,^4. In 
principle, we solve the four conditions for the four variables Oq, Oq, £1, £2, treating 
the symmetry breaking coupling constants £3, £4, ri, . . . ,r4 as known. 

The constant £3 is varied in the range specified in ( [7.81 ). Concerning £4, we rely 
on the result for the scalar radius given in ( [7.7D , thus in effect replacing the input 
variable £4 by {r'^)s- The analysis then involves a fifth condition: the relation 
(|7.5| ), which expresses the scalar radius in terms of effective coupling constants. 

If all of the input variables are taken at their central values, the representation 
for the moments given in appendix ^ can be used. The solution of the resulting 
system of numerical equations occurs at the values quoted in table |l], first row. 
The next four rows indicate the sensitivity to the input used for (r^)^, £3, to 
the uncertainties in the symmetry breaking coupling constants r„ of 0{p^), and 
to those in the experimental information used when solving the Roy equations. 
The details of the error analysis that underlies these numbers are described in 
appendix 0. 

The table shows that the uncertainties in the prediction for Oq and are 
dominated by those from £3. In particular, the result for the S'-wave scattering 
lengths is not sensitive to the contributions from the coupling constants occurring 



19 



at Oijf'). Adding up the uncertainties due to these and to the experimental input 
in the Roy equations, we arrive at 

al= 0.220 ±0.001 + 0.027 A^2 -0.0017 A4, (11.1) 
al = -0.0444 ± 0.0003 - 0.004 A^2 - 0.0004 A^s , 

where A ^2 and Ai^ are defined by 

(r^)^ = 0.61 irn^il + A,2) , £3 = 2.9 + A^ . 

Our final result for the S'-wave scattering lengths follows from this representation 
with the estimates for Ar2, Ai^ given in ( |7.8| ), ([7.7|), and reads 

a° = 0.220 ±0.005, = -0.0444 ± 0.0010 , (11.2) 

2a° - 5al = 0.663 ± 0.006 , 0° - = 0.265 ± 0.004 . 

Expressed in terms of the coefficients Cq, Ci, C2, this result corresponds to 

Co = 1-096 ± 0.021 , Ci = 1.104 ±0.009, C2 = 1.115 ± 0.022 . (11.3) 



12 Discussion 

The terms omitted in the chiral perturbation series represent an inherent limita- 
tion of our calculation. The matching must be done in such a manner that these 
are small. In contrast to a matching at threshold - that is, to the straightforward 
expansion of the scattering lengths - our method fulfills this criterion remarkably 
well: We are using the expansion in powers of the quark masses only for the 
coefficients Cq, Ci and C2, while the curvature generated by the unitarity cut is 
evaluated phenomenologically. As discussed in section |], the infrared singulari- 
ties occurring in the expansion of these quantities have remarkably small residues. 
Indeed, truncating the expansion of C„ at order 1, m and m^, respectively and 
solving equation ( |6.5| ) in the corresponding approximation, we obtain 

= 0.197 ^ 0.2195 ^ 0.220, 
al = -0.0402 ^ -0.0446 -0.0444 , (12.1) 
2a^Q-5al= 0.594 ^ 0.662 0.663, 

indicating that the series converges very rapidly. For this reason, we expect the 
contributions from yet higher orders to be entirely negligible. 

The rapid convergence of the series is a virtue of the specific method used 
to match the chiral and phenomenological representations. To demonstrate this, 
we briefiy discuss the alternative approach used in refs. 0, |^ , where the results 
for the various scattering lengths and effective ranges are obtained by directly 
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evaluating the chiral representation of the scattering amphtude at threshold. 
Keeping the values of the effective coupling constants fixed at the central values 
and truncating the series at order m, and m^, we obtain the sequence 



2 - 5 



0.159 ^ 0.200 ^ 0.216, 
-0.0454 -0.0445 -0.0445 
0.545 0.624 0.654. 



;i2.2) 



The first terms on the right correspond to Weinberg's formulae. The second and 
third terms are in agreement with the old one loop results of ref. |^ and the 
two loop results of ref. 22, 44 1, respectively. As indicated by the difference 
between the second and third terms, the corrections of 0{p^) are by no means 
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Figure 1: Constraints imposed on the S'-wave scattering lengths by chiral sym- 
metry. The three full circles illustrate the convergence of the chiral perturbation 
series at threshold, according to eq. (12.2). The one at the left corresponds to 
Weinberg's leading order formulae. The error ellipse represents our final result. 
The other elements of the figure are specified in the text. 

negligible for a matching at threshold. This is illustrated in fig. |I|, where the three 
full circles correspond to the sequence (|12.2|) . The triangle at the right and the 



shaded rectangle indicate the central values and the uncertainties quoted in the 



1979 compilation of ref. [^. The triangle and the diamond near the center of the 



figure correspond to set I and set II of ref. ||^, respectively. The elhpse represents 
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the 68% confidence contour of our final result in eq. ( |11.2| ). The details of the 
error analysis that underlies this result are described in appendix ^ 

The reason why the straightforward expansion of the scattering lengths in 
powers of the quark masses converges rather slowly is that these represent the 
values of the amphtude at threshold, that is at the place where the branch cut 
required by unitarity starts. The truncated chiral representation does not de- 
scribe that singularity well enough, particularly at one loop, where the relevant 
imaginary parts stem from the tree level approximation. 

If the effective coupling constants are the same, the only difference between 
our method and a matching at threshold is the one between the functions W\s) 
and U^{s). In particular, the results for aj], Oq only differ because the numerical 
values of W^{s) and U^{s) at s = AM^ are not the same. As mentioned above, 
the difference between the two sets of functions affects the scattering amplitude 
only at 0{p^) and beyond. Numerically, however, it is not irrelevant which one of 
the two is used to describe the effects generated by the unitarity cuts: While the 
functions W\s) account for the imaginary parts of the S- and P- waves to the 
accuracy to which these are known, the quantities U\s) represent a compara- 
tively crude approximation, obtained by evaluating the imaginary parts with the 
one-loop representation. 

13 Correlation between and 

As mentioned earlier, the main difference between Generalized Chiral Pertur- 
bation Theory and the standard one used in the present paper resides in the 
coupling constant £3. Apart from that, the formulae are identical - only the 
bookkeeping for the chiral power of the quark mass matrix is different. In par- 
ticular, the relation between the scalar radius and the coupling constant £4 also 
holds in that framework, but there is no prediction for the S'-wave scattering 
lengths Og and Oq, because these involve the coupling constant £3. The fact that 
£4 is strongly constrained by the value of the scalar radius implies, however, that 
there is a strong correlation between Oq and Oq, independently of whether the 
quark condensate is the leading order parameter: Apart from higher order cor- 
rections, both of these are controlled by the same parameter £3. The dependence 
is approximately described by the parabolae 

= 0.225 - 1.6 • 10"^ £3 - 1.3 ■ 10'^ (Isf , (13.1) 
al = -0.0433 - 3.6 ■ 10"^ I3 - 4.3 ■ 10"^ . 

which are displayed in fig. Note that the interval shown far exceeds the range 

^If £3 is large, the symmetry breaking effects generated by tire quark masses are larger than 
in the standard framework, so that a reordering of the series that gives these more weight is 
called for. 
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Figure 2: S'-wave scattering lengths as functions of 



relevant for the standard picture, which is indicated by the vertical bar. 
Eliminating the parameter £3, we obtain a correlation between Oq and a^: 



-.0444 ± .0008 + .236 (a^ - .22) 
-.61 (a° - .22)2 - 9.9 (aj^ - .22)^ 



;i3.2) 



The error given accounts for the various sources of uncertainty in our input - 
evaluating these as described in appendix we find that they are nearly in- 
dependent of Qq. The correlation is indicated in fig. ^ The values of Og and 
Oq are constrained to the region between the two dash-dotted lines that touch 
the error ellipse associated with the standard picture. As discussed in ref. 0, a 
qualitatively similar correlation also results from the Olsson sum rule ||^ - the 
two conditions are perfectly compatible, but the one above is considerably more 
stringent. Fig. |I]also shows that for Oq < 0.18, or £3 > 25, the center of the region 
allowed by the correlation falls outside the universal band, which is indicated by 
the tilted lines. The same happens on the opposite side, for Oq > 0.28, £3 < —54. 
Since the Roy equations only admit solutions if the two subtraction constants Og 
and Oq are in the universal band, exceedingly large values of £3 are thus excluded. 
Note also that the correlation implies an upper bound on the 1 = 2 scattering 
length: < -0.035. 

The correlation between and can be used, for instance, to analyze the 
information about the phase difference 5q — 6l obtained from the decay K — »• 
vrvre u. At the low energies occurring there, this difference is dominated by the 
contribution oc Og from the / = S'-wave scattering length. The relation (|13.2| ) 
allows us to correct for the higher order terms of the threshold expansion: The 
phase difference can be expressed in terms of the energy and the value of Og, up 
to very small uncertainties. This is illustrated in fig. S The center of the three 
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Figure 3: Phase relevant for the decay K ttttcz/. The three bands correspond 
to the three indicated values of the S-wave scattering length a^. The uncertainties 
are dominated by those from the experimental input used in the Roy equations. 
The triangles are the data points of Rosselet et al. 
represent the preliminary E865 results . 



28|, while the full circles 



narrow bands shown is obtained by fixing the value of Oq with the correlation 
( |13.2|) and inserting the result in the numerical parametrization of the phase 
shifts in appendix D of ref. 0. At a given value of Oq, the uncertainties in the 
result for the phase difference 5q{s) — 5\{s) are dominated by the one in the 
experimental input used for the / = 5'-wave. Near threshold, the uncertainties 
are proportional to (s — 4M^)^''^ - in the range shown, they amount to less than 
a third of a degree. While the data of Rosselet et al. ||2^ are consistent with all 
three of the indicated values of Oq, the preliminary results of the E865 experiment 
^ are not. Instead they beautifully confirm the prediction 

0.218, with = 5.7 for 5 



at Brookhaven 



The best fit to these data is obtained for Oq 



degrees of freedom. As pointed out in ref. |^ , the correlation ( p.3.2| ) can be used 
to convert data on the phase difference into data on the scattering lengths. For 

we refer to [H6 



a detailed discussion of the consequences for the value of Og, 



47]. 



14 Results for £i and £2 

The effective coupling constants of £4 enter the chiral perturbation theory repre- 
sentation of the scattering amplitude and of the scalar form factor only as correc- 
tions, so that our results for these are subject to significantly larger uncertainties 
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than those for Oq, Oq. According to table |1], we obtain 

£1 = -0.4 ±0.6, £2 = 4.3 ±0.1. (14.1) 

The noise in the symmetry breaking couphngs r„ of and the one in the Roy 
equation input yield comparable contributions, while those from the other entries 
are negligibly small. The corresponding error ellipse is shown in fig. ^. 
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Figure 4: Values of the coupling constants ^l and ^2- The shaded ellipse shows 
the result of our calculation. The rectangles indicate the ranges quoted in refs. ||^, 
^ and |2^. The triangle and the diamond correspond to set I and set II of 



respectively. The cross represents the resonance saturation estimate of ref. [Q. 
The full circle is the result obtained by matching at one loop and the thin ellipse 
close to it represents the uncertainties in the effective one loop couplings , ■ 



In order to investigate the uncertainties due to the neglected higher order 
terms, we again compare this with what is found if the phenomenological rep- 
resentation is matched to the one loop approximation of the chiral perturbation 
series. For the central values of the input parameters, the solution of the match- 
ing conditions then occurs at £1 = —1.8, I.2 = 5.4 : The two loop effects shift the 
one loop result by about ±1.4 and —1.1 units, respectively. The shift arises from 
the fact that the expansion of the coefficients C3 and C4 contains very strong in- 
frared singularities at first nonleading order. Analogous contributions also occur 
in ci and C2, at next-to-next-to-leading order, but in the combinations Co, Ci, 
C2 that matter for the determination of the scattering lengths, these singularities 
only generate very small effects: In these quantities, the contributions of order 
amount to less than 1%. We conclude that, unlike the result for Oq, Oq, where the 
uncertainties from the neglected higher order terms are tiny, the one for ii and 
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^2 is sensitive to these. Although we expect the corresponding contributions to 
be small compared to the first order shift given above, they might be of the same 
order as those from the uncertainties in our input - we do not offer a quantitative 
guess. 

The couplings ii and £2 are quark mass independent, whereas the physical 
quantities used to estimate their values incorporate quark mass effects. As a 
result of this, it is problematic to rely on phenomenological determinations based 
on the one loop approximation when analyzing quantities at two loop order. 
The large infrared singularities that accompany the contributions from ii and £2 
are automatically accounted for in the two loop representation, but are missing 
in the framework of a one loop calculation - in the phenomenological analysis, 
their contributions are lumped into those from the coupling constants. As an 
illustration, we mention the set I of couplings introduced in 0], that uses the 
one- loop values for ii and £2, but leads to D-wave scattering lengths that do not 
agree well with the values extracted from experiment, as was first pointed out in 
ref. For a detailed discussion of this issue, we refer to 



We now show that, once the shift in the values of ii, £2 is accounted for, the 
one and two loop representations for the coefficients Ci, . . . ,04 become nearly the 
same, so that the results obtained by matching the phenomenological represen- 
tation with the chiral one at two loop level nearly coincide with those found in 
the one loop approximation. The infrared singularities responsible for that shift 
are those contained in the coefficients 63, 64. If we solve the expressions for these 
coefficients in one loop approximation, we obtain 

£f ^3ih-h) + ^, £f^6h + ^. (14.2) 

The expansion of these quantities in powers of the quark masses starts with 
£f = In + 0{^). The infrared singularities generated by the two loop graphs 
show up in the terms of order ^, in particular through contributions proportional 
to P = In^ (/iVM2), which are very important numerically. Accounting for the 
uncertainties in our input, we obtain 

if = -1.9 ± 0.2 , £f = 5.25 ± 0.04 . (14.3) 

The comparison with the values £1 = —1.8, £2 = 5.4, found when matching 
at one loop, shows that the couplings relevant in the context of the one loop 
approximation may indeed be characterized in this manner (compare fig. ^, where 
the values for £1, £2 obtained at one loop are indicated by the full circle, while 
the thin ellipse corresponds to the above numerical result for £f, £f). 

Now comes the point we wish to make: We may also evaluate the one loop 
formulae ( |6.1|) for ci, C2, replacing £1, £2 by the above effective values. The 
outcome differs from what is obtained with the two loop formulae only by a 
fraction of a percent - the difference is in the noise of the two loop result. In this 
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sense, the main effect of the infrared singularities in the two loop graphs amounts 
to a shift in the values of the coupling constants Zi, £2- This explains why the 
matching conditions used in the present paper yield very accurate results for the 
5'-wave scattering lengths already at one loop, while the corresponding results 
for these two couplings are off. 

The literature contains quite a few determinations of the coupling constants 
£1 and £2 that are based on the one loop approximation of chiral perturbation 



theory - |]T6|, starting with the estimates £1 = —2.3 ±3.7, £2 = 6.0 ±1.3 given 
in ref. 0, which are perfectly consistent with our result for £f-^, €2 ■ Note that, in 
the case of £2-, the shift generated by the two loop graphs takes the result outside 
the quoted range (as stated in ref. 0, that range only measures the accuracy 
to which the first order corrections can be calculated and does not include an 
estimate of contributions due to higher order terms). 

The results for the effective coupling constants obtained by Girlanda et al. 



read £x = -0.37±0.95± 1.71, £2 = 4.17±0.19±0.43. The first error comes from 
the evaluation of the integrals over the imaginary parts, while the second reflects 
the uncertainties in the contributions from the couplings of Cq. Our results in 
eq. ( |14.1 ) confirm these numbers, with substantially smaller errors - we repeat, 
however, that these only account for the noise seen in our calculation. 

Amoros, Bijnens and Talavera |^2|] have extracted values for the coupling 
constants of £4 from their two loop analysis of the K^^ form factors - which is 
based on SU(3)flXSU(3)£, chiral perturbation theory - and obtain £1 = 0.4 ±2.4, 
£2 = 4.9 ±1.0. Fig. ^ shows that these are perfectly consistent with ours. As these 
authors are relying on the one loop relations between the coupling constants L„ 
of that framework and the couplings £n relevant for SU(2)j5xSU(2)^, the results 
are accompanied by comparatively large errors. 



15 Values of £4, and tq 

For the central values of the input, the matching conditions lead to £4 = 4.39 
(first row in table |l]). The uncertainties in this number due to the various sources 
of error are dominated by the one in the scalar radius and the noise in the 
symmetry breaking coupling constants ri, r2, r^, r4, of Cq. In order to 
estimate the uncertainties due to the higher order effects that our calculation 
neglects, we compare the above two loop result with the value £4 = 4.60, obtained 
by truncating the chiral representation for the scalar radius at leading order. The 
comparison shows that the shift generated by the two loop contributions is of the 
same size as the one due to the uncertainty in the scalar radius. Those from yet 
higher orders are expected to be significantly smaller, so that the uncertainty in 
the final result is dominated by the sources of error listed in the table. The net 
result reads 

£4 = 4.4 ±0.2. (15.1) 
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The number is consistent with the one loop estimate £4 = 4.3 ± 0.9, given in 
ref. 0. The infrared singularities that accompany the coupling constant £4 are 
much weaker than those occurring together with ii, £2- The same is true also for 
£3, where the uncertainties are much too large for such effects to matter at all. 
The above result confirms the value £4 = 4.4 ± 0.3, obtained by Bijnens, 



Colangelo and Talavera from a comparison of the two loop representation 
with the dispersive result of the scalar radius, but this was to be expected, because 
the input used in the two evaluations is nearly the same. 



In the framework of the calculation mentioned in section 14, Amoros, Bijnens 



and Talavera |^ obtain £4 = 4.2 ± 0.18, also consistent with our result (as 
emphasized by these authors, the error bar does not account for the uncertainties 
due to higher order effects, which in their approach are quite substantial). 

The coupling constants Vn = (4vr)'' r^(/i) are scale dependent. We could intro- 
duce corresponding scale independent quantities, analogous to the terms In used 
for the coupling constants of £4. The scale dependence is rather complicated, 
however, because it is quadratic in In/i. We instead quote the values obtained 
for fi = Mp = 0.77 GeV. Our analysis does not shed any light on the symmetry 
breaking coupling constants ri, . . . , r4, which belong to the input of our calcula- 
tion, but we can determine and rg, from the matching conditions for C5 and Cq 

- we did not yet make use of these. Numerically, we find: 

f5 = 3.8±1.0, f6 = l-0±0.1. (15.2) 

Table shows that the noise seen in our calculation is dominated by the one in the 
estimates for the symmetry breaking coupling constants ri, . . . , r4. Note that the 
error bars do not account for the uncertainties due to higher order contributions 

- our evaluation does not give us any handle on these. 

The resonance estimates of refs. |4^ offer a test: They lead to 



f5^2.7, ^6^0.75, (15.3) 
and thus corroborate the outcome of our analysis, both in sign and in magnitude. 



In fact, as pointed out by Ecker [48], the estimates 



-0.7, £2^5.0, 4^1-9, £4^3.7, (15.4) 



obtained from resonance saturation of sum rules |3^, are perfectly consistent 
with the numbers found at two loop accuracy. We conclude that there is good 
evidence for the picture drawn in ref. to be valid: The values of all of the ef- 
fective coupling constants encountered in the two loop representation of the scat- 
tering amplitude are consistent with the assumption that these are dominated 
by the contributions from the singularities due to the exchange of the lightest 
non-Goldstone states. Admittedly, this assumption does not lead to very sharp 
values, because the separation of the resonance contributions from the continuum 
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underneath is not unique. The problem manifests itself in the scale dependence 
of the coupling constants - resonance saturation can literally hold only at one 
particular scale. Also, it is not a straightforward matter to formulate the reso- 
nance saturation hypothesis for singularities due to the exchange of particles of 
spin two or higher |19|, Even so, we consider it important that the values 
found for the coupling constants are within the noise inherent in the assumption 
that, once the poles and cuts due to the Goldstone bosons are removed, the low 
energy behaviour of the scattering amplitude is dominated by the singularities 
due to the remaining states. Since these remain massive in the chiral limit, their 
contributions to the chiral expansion are suppressed by powers of momenta or 
quark masses, but they do show up at nonleading orders. 



16 The coefficients 61, . . . , 



The matching conditions ( |4.2| ) express the coefficients Cn of the chiral represen- 
tation in terms of the S'-wave scattering lengths and moments of the imaginary 
parts. Inserting the numerical representation for the dependence of the moments 



on the scattering lengths and comparing the result with eq. (|B.2|) , we obtain the 



following representation for the coefficients introduced in ref. ||^: 

61 = -.1 ± .1 - 21 Aa° + 1670 Aal + 9 (Aa°)2 + 96 Aa°Aa^ - 972 (Aalf, 

62 = 8.2 ± .4 + 179 Aa° - 602 Aal - 135 (Aa")^ + 315 Aa^Aa^ - 65 (Aalf, 

63 = -.41 ± .06 + 3.5 Aal - 12.9 Aal + 7 (^«o)^ - 30 Aa^Aa^ + 40 {Aalf, 

64 = .71 ± .01 + 1.3 Aal - 4-1 Aag - {^alf - 4 Aa°Aag + 25 (Aalf, (16.1) 

65 = 2.99 ± .35 + 32.6Aa° - 97.0 Aal + 104(Aa°)2 - 451Aa°Aa^ + 602{Aaiy, 
Iq = 2.18 ± .01 + 7.2 Aal - 28.4 Aal - 3 (Aflo)^ + 9 AalAal - 62 (Aalf, 

with Aal = '^0 ^ 0.225, Aal = '^0 + 0.03706. The error bars indicate the uncer- 
tainties in the outcome due to those in the experimental input used when solving 
the Roy equations. The representation holds for arbitrary values of the scattering 
lengths in the vicinity of the point of reference. Inserting our results from ( 11. 2| ) 
and adding errors quadratically, we finally obtain 

61 = -12.4 ± 1.6, 62 = 11.8 ±0.6, 63 = -0.33 ±0.07, 

64 = 0.74 ±0.01, 65 = 3.58 ±0.37, fog = 2.35 ± 0.02 . (16.2) 

We emphasize that the error bars only indicate the noise seen in our evaluation. 
In 61, . . . ,64, the two loop representation does account for the contributions of 
next-to-leading order, but in the case of 65, Bq, it only yields the leading terms 
- these quantities are particularly sensitive to the neglected higher orders. 

The above results may be compared with the values found in the literature. 
Girlanda, Knecht, Moussallam and Stern [^] work within generalized chiral per- 
turbation theory and do not have a prediction for the magnitude of the coefficients 
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Figure 5: Result for 63 and 64. The errors in our result are dominated by those 
in the experimental input used when solving the Roy equations: The nearly 
degenerate ellipse indicates the result obtained if these could be ignored. The 
rectangles correspond to the values quoted in refs. 
marks the one obtained in ref. set II. 



181 and 24 , while the diamond 



bi and 62, because the corresponding expressions contain the two free parame- 
ters a and /3. In their framework, the analogs of the constants 63, . . . , fee are 
denoted by Ai, . . . , A4. The explicit relation between the two sets of quantities is 
given in eq. ([A.l|) . In our notation, the numerical values of ref. correspond 
to 63 = -0.56 ± 0.37, 64 = 0.737 ± 0.039, 65 = 3.25 ± 1.50, h = 2.42 ± 0.22 
and are perfectly consistent with our results, where the errors are smaller. The 
result for bi and 62, obtained above within the standard framework, amounts to 
a prediction for the magnitude of a and (3. Numerically, we obtain 



a 



1.08 ±0.07. 



/3 = 1.12 ±0.01 



(16.3) 



Wanders ITS 



has obtained values for the coefficients 63, 64 and feg from mani- 
festly crossing symmetric dispersion relations. Matching the chiral and dispersive 
representations at the center of the Mandelstam triangle, he obtains the values 
63 = -0.403 ± 0.032, h = 0.719 ± 0.024, be = 2.29 ± 0.075, which are also 
consistent with our numbers. Note that the quoted errors only account for the 
uncertainties arising from the procedure used in ref. [0 and do not cover those 
in the input. Fig. ^ shows that, in the case of 63, the experimental input in the 
Roy equations represents the dominating source of error. 



Amoros, Bijnens and Talavera |^2[ have determined the coefficients 6„ on the 
basis of their analysis of the K^^ form factors, referred to earlier. The results for 
the coefficients &3, . . . ,65 are accompanied by rather large errors and we do not 
list these here, but merely note that the central values in eq. ( |16.2| ) are within the 
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Figure 6: Result for 65 and 65 • The strip between the two horizontal lines corre- 



sponds to the value for of Wanders |18 



quoted range, in all cases. For the first two terms, however, Amoros et al. arrive 
at comparatively accurate values, bi = —10.8 ± 3.3, 62 = 10.8 ± 3.2, which are 
also perfectly consistent with those in eq. ( |16.2| ). The fact that, in their analysis, 
the remaining coefficients are subject to large uncertainties, also manifests itself 
in column C of table 0: The error bars in the first five rows of the table, Qq ■ ■ ■ a\, 
are much smaller than those in the remainder. 



17 S— and P— wave phase shifts 

For the reasons discussed in detail in ref. the two S-wave scattering lengths 
are the essential parameters in the low energy domain. The result in eq. ( |11.2| ) 



specifies these to within very small uncertainties. In particular, we can now work 
out the phase shifts of the S- and P-waves on this basis, using the Roy equation 
analysis of 0. The available experimental information for the imaginary parts 
above ^/so = 0.8 GeV, as well as the scattering lengths a^, are used as an 
input, while the output of the calculation consists of the phases Sq{s), Sl{s) 
and Sq{s), in the region below Sq. In view of the two subtractions occurring in 
the Roy equations, the behaviour of the imaginary parts above 1 GeV has very 
little infiuence on the behaviour of the solutions below 0.8 GeV. Also, there is 
a consistency check: In the region above Sq, the output must agree with the 
input. For the values of the scattering lengths required by chiral symmetry, this 
condition is indeed met. In fact, the solution of the Roy equations closely follow 
the input, within the rather broad range of variations allowed for the imaginary 
parts in ref. 0. This also means that the Roy equations do not strongly constrain 
the behaviour of the phases above 0.8 GeV. 

The result is shown in figs. |^, ^ and ^. For comparison, these figures also 



31 




0.3 0.4 0.5 0.6 0.7 0.8 

E(GeV) 



Figure 7: 1 = S-wave phase shift. The full line results with the central values 
of the scattering lengths and of the experimental input used in the Roy equations. 
The shaded region corresponds to the uncertainties of the result. The dotted lines 
indicate the boundaries of the region allowed if the constraints imposed by chiral 
symmetry are ignored |^ . The data points are from refs. and ||5^ . 



show the data points of the phase shift analyses given by Hyams et al. 
Protopopescu et al. the solutions A and B of Hoogland et al. (ACM) 



and the one of Losty et al. as well as the P-wave phase extracted from the 
data on the reactions e'^e~ —>■ tt+tt" and r —>■ wkit. For further information on 
the S-wave phase shifts, we refer the reader to l5F 
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The three central curves are described by the parametrization [p5 



(17.1) 
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Figure 8: P-wave phase shift. The phase of the pion form factor is also shown, 
but it can barely be distinguished from the central result of our analysis. The 



data points are from refs. [^] and 
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17.2) 



in units of M^. In particular, the constants A| represent the scattering lengths 
of the three partial waves under consideration, while the are related to the 
effective ranges. 

The parameters specify the value of s where passes through 90°: 



36.77 Mi 



30.72 m: 



-21.62 Mi 



(17.3) 



In the channels with / = 0,1, the corresponding energies are 846 MeV and 
774 MeV, respectively (the negative sign of Sq indicates that in the 1 = 2 channel, 
which is exotic, the phase remains below 90°). 
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Figure 9: 1 = 2 S-wave phase shift. The full line results with the central 
values of the scattering lengths and of the experimental input used in the Roy 
equations. The shaded region corresponds to the uncertainties of the result. The 
data points represent the two phase shift representations of the Aachen-Cern- 
Munich collaboration ||51| and the one of Losty et al. 152 1 



The value of the phase difference 6q — 6q a.t s = M]^ is of special interest, in 
connection with the decays K tttt. In particular, the phase of e'/e is determined 
by that phase difference. Our representation of the scattering amplitude allows 
us to pin this quantity down at the 3% level of acuracy: 

50(M|.o) - 5o'(^io) = 47.7° ± 1.5° . (17.4) 

We add two remarks concerning the comparison with the P-wave phase shift 
extracted from the e~^e~ and r data. First, we note that the agreement at 0.8 GeV 
is enforced by our approach: In the Roy equation analysis, the value of the phase 
shift at that energy represents an input parameter and we have made use of those 
data to pin it down. Once that is done, however, the behaviour of the phase shift 
at lower energies is unambiguously fixed: Chiral symmetry determines the two 
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subtraction constants, so that the solution of the Roy equations becomes unique. 
In other words - disregarding the small effects due to the uncertainties in the 
input of our analysis, which are shown in fig. ^ - there is only one interpolation 
between threshold and 0.8 GeV that is consistent with the constraints imposed 
by analyticity, unitarity and chiral symmetry. Figure ^ shows that the predicted 
curve indeed very closely follows the phase extracted from the e'^e~ and r data. 
This confirms the conclusions reached in ref. ||56[| . 

Actually, the figure conceals a discrepancy in the threshold region, where the 
phase is too small for the effect to be seen by eye: Evaluating the P-wave scatter- 
ing length with the Gounaris-Sakurai parametrization of the form factor given in 
ref. ||5^ (the curve shown in the figure), we obtain a result that is smaller than the 
value for a\ in table 0, by about 10%, that is by many standard deviations of our 
prediction. The discrepancy is in the noise of the data on the form factor: There 
is little experimental information in the threshold region, so that the behaviour 
of the form factor is not strongly constrained there. Indeed, there are alternative 
parametrizations that also fit the data, but have a distinctly different behaviour 
near and below threshold. Even parametrizations with unphysical singularities at 
s = 0, such as the one proposed in provide decent fits in the experimentally 
accessible region. In this respect, the present work does add significant informa- 
tion about the P-wave phase shift, as it predicts the behaviour near threshold, 
within very narrow limits. 

18 Poles on the second sheet 

The partial wave amplitude t\{s) contains a pole on the second sheet. Denoting 
the pole position by s = (Mp — | Ep)^, we obtain 

Mp = 762.4 ± 1.8 MeV, Ep = 145.2 ± 2.8 MeV . (18.1) 

Note that the values quoted for the "mass" often represent the energy where 
the real part of the amplitude vanishes - in contrast to the position of the pole, 
that value is not independent of the process considered. As the scattering is 
approximately elastic there, the corresponding mass is the energy where the phase 
shift goes through 90°. Eor the P-wave, this happens at 

rup = 773.5 ±2.5 MeV. 

The real part of the pole position is smaller than the energy where the phase 
shift passes through 90°, by about 10 MeV. The uncertainty in Ep is significantly 
smaller than the error bar quoted in 0: The constraints imposed on the scattering 
amplitude by the low energy theorems of chiral symmetry also allow a better 
determination of the width. 

The 1 = S-wave also contains a pole on the second sheet. The uncertainties 
in the pole position are considerably larger than in the case of the p, because the 
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singularity is far from the real axis. Also, the uncertainties in the phase shift 
are somewhat larger here. Varying the input parameters as well as the analytic 
form of the representation used for tg, we find that the pole occurs in the region 
^/s = (470 ±30) - i (295 ± 20) MeV, while the phase passes through 90° at 
v/i = 844 ± 13 MeV. 

There is no harm in calling this an unusually broad resonance, but that sheds 
little light on the low energy structure of the scattering amplitude. In particular, 
it should not come as a surprise if the values for the mass and width of the 
resonance, obtained on the basis of the assumption that the pole represents the 
most important feature in this channel, are very different from the real and 
imaginary parts of the energy at which the amplitude actually has a pole - there 
is more to the physics of the S'-wave than the occurrence of a pole far from the 



real axis. A collection of numbers concerning the pole position is given in |^ 
and for a recent review of the abundant literature on the subject, we refer to 
A recent discussion in the framework of the NN interaction is given in [HT|. 



We add a remark concerning the physics behind the pole in - admittedly, 
the reasoning is of qualitative nature. In the chiral limit, current algebra predicts 
to = s/IGtcF^: The amplitude vanishes at threshold, but the real part grows 
quadratically with the energy, so that the imaginary part rises with the fourth 
power. The rapid growth signals the occurrence of a strong final state interaction. 
In order to estimate the strength of the corresponding branch cut, we invoke the 
inverse amplitude method, replacing the above formula by tg = s/ (IGttF^ — is). 
The virtue of this operation is that, while it retains the algebraic accuracy of 
the current algebra approximation, it yields an expression that does obey elastic 
unitarity. The formula shows that, in this approximation, the amplitude contains 
a pole at ^/s = 16 n = 463 — 2 463 MeV, indeed not far away from the 
place where the full amphtude has a pole. 

The physics of the P-wave is very different, because the unitarity cut gen- 
erated by low energy vrvr intermediate states is very weak. Repeating the above 
exercise for t\, one again finds a pole with equal real and imaginary parts, but it 
is entirely fictitious, as it occurs at 1.1 — i 1.1 GeV, far beyond the region where 
current algebra provides a meaningful approximation. The occurrence of a pole 
near the real axis cannot be understood on the basis of chiral symmetry and 
unitarity alone. 

In the framework of the effective theory, the difference manifests itself as 
follows. While the unitarity corrections account perfectly well for the low energy 
behaviour of the imaginary parts, the presence of the p only shows up in the 
values of the effective coupling constants £1,^2 and i^. There is no such pole in 
t\, for instance, if the underlying theory is identified with the linear a-model, 
and the values of those coupling constants are then very different |0. In this 
sense, the pole in t\ reflects a special property of QCD, while the one in tg can 
be understood on the basis of the fact that chiral symmetry predicts a strong 
unitarity cut: The pole position is related to the magnitude of F^. 
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19 Threshold parameters 





A 


B 


C 


D 


E 


units 


hi 


.220 ± .005 
.276 ± .006 


.216 
.268 


.220 ± .005 
.280 ± .011 


.24 ± .06 
.26± .02 


.26 ± .05 
.25± .03 




al 

u 

hi 


-.444 ± .010 
-.803 ± .012 


-.445 
-.808 


-.423 ±.010 
-.762 ±.021 


-.36± .13 
-.79± .05 


-.28± .12 
-.82± .08 


10-1 

lO-iM-2 


a\ 
h\ 


.379 ± .005 
.567 ± .013 


.380 
.537 


.380 ± .021 
.58 ±.12 


.37± .02 
.54± .04 


.38± .02 


lO-iM-2 
lO-^M-^ 


al 
hi 


.175 ± .003 
-.355 ± .014 


.176 
-.343 


.22 ±.04 
-.32 ± .10 


.17± .01 
-.35± .06 


.17± .03 


10-2M-4 


al 
hi 


.170 ± .013 
-.326 ± .012 


.172 
-.339 


.29 ±.10 
-.36 ±.9 


.18± .08 
-.34± .07 


.13± .3 


10-^M-^ 
10-3M-6 


a\ 
h\ 


.560 ± .019 
-.402 ± .018 


.545 
-.312 


.61 ± .11 
-.36 ±.02 


.58± .12 
-.44± .14 


.6± .2 


lO-^M-6 
lO-^M-8 



Table 2: Threshold parameters. Our results are listed in column A. The numbers 
in the next two columns are obtained by evaluating the chiral representation at 
threshold: The entries under B follow from our values of the effective coupling 
constants, while those under C are taken from ref. E^. Column D gives the 
outcome of a Roy equation analysis that does not invoke chiral symmetry ||^, 
while E contains the old "experimental" values 



The scattering lengths of the partial waves with £ > 1, as well as the effective 
ranges (also of those of the ^-waves) can be expressed in terms of sum rules over 
the imaginary parts p2[. The corresponding numerical values are listed in the 
table 0, together with the S'-wave scattering lengths. Column A indicates our 
final results, obtained by matching the phenomenological and chiral represen- 
tations in the subthreshold region and using the Roy equations to evaluate the 
amplitude and its derivatives at threshold. In column B, we quote the numbers 
obtained from a direct evaluation of the two loop representation at threshold, 
using our central values for the effective coupling constants - this amounts to 
truncating the expansion of the threshold parameters in powers of and m^. 
Column C lists the results of ref. 



n 



where the amplitude is also expanded at 
threshold, but the coupling constants are determined on the basis of a two loop 
analysis of the Kei form factors (see the next section for a comment concerning 
these entries). The comparison of the columns A, B and C clearly shows that two 
loop chiral perturbation theory works very well in describing both tttt scattering 
and Kf,i decays. For the reasons given in section ^ the method described in 
the present paper yields the smallest error bars. In fact, it is quite remarkable 
that the results for the effective range 63 in columns B and C represent a decent 
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estimate: In this case, only the infrared singularities occurring in the expansion 
in powers of the quark masses contribute. For comparison, column D lists the 
values of ref. 0, which are obtained by analyzing the available data with the 
Roy equations and do not invoke chiral perturbation theory. Finally, column E 
contains the values of the compilation in refs. . 



20 Quark mass dependence of and F^^ 

The dependence of physical quantities on the quark masses is of interest, in 
particular, for the following reason: By now, dynamical quarks with a mass of 
the order of the physical value of are within reach on the lattice, but it is 
notoriously difficult to equip the two hghtest quarks with their proper masses. 
Invariably, the numerical results obtained for the physical values of m„ and 
rely on an extrapolation of numerical data. For a recent, comprehensive review 
of lattice work on the light quark masses, we refer to ||63l . 



In this connection, chiral perturbation theory may turn out to be useful, 
because it predicts the mass dependence in terms of a few constants: the couphng 
constants of the effective Lagrangian. Above, we have determined some of these 
and we now discuss the consequences for the dependence of M^, F.„,aQ, on the 
masses of the two lightest quarks: We keep rus fixed at the physical value and 
set rriu = = m, but vary the value of m in the range < m < | (at the 
upper end of that range, the pion mass is about 500 MeV). 

In the preceding sections, we have expressed all of the quantities in terms 
of the physical pion mass and the physical decay constant, using the ratio 
as an expansion parameter. Also, the logarithmic infrared singularities were 
normalized at the scale M^r. In particular, the coupling constants In contain a 
chiral logarithm with unit coefficient, so that they may be represented as 

C-ln^, (20.1) 

where A„ is the intrinsic scale of in and is independent both of m and of the 
running scale yU. In order to explicitly exhibit the quark mass dependence, we 
replace by the variable = 2mB and also normalize the chiral logarithms 
at the scale M, trading the quantities ^ and In for 

M2 . A 2 

x = 4 = ln^, (20.2) 



respectively. According to (|7. 1|) , (|7.2| ) the two sets of variables are related by 

x = a^+iai3+ 4^4) + oien , 4 = ^ - k ^3 + o^e) (20.3) 

e = X{1 - i X (4 + 4^4) + 0(x')} , In = L + kx is + 0{X^) . 
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The expansions of and in powers of m are known to next-to-next-to- 
leading order [^4|, ^. In the above notation, the exphcit expressions may be 



written in the form |66|: 



= {1 - i a; 4 + f x^i^, + x^h, + Oix^) , 
= F {1 + X £4 - I + x^kp + O(x^) , 

4 = ^ (28 li + 32 4 -94 + 49) , (20.4) 
51 

1 = ^ (14^\ + 164 + 64 - 64 + 23) , 

In this representation, the infrared singularities are hidden in the scale invariant 
quantities £1, . . . ,£4. Those generated by the two loop graphs have been com- 
pleted to a square. The normalization of the auxiliary quantities 4/, 4 is chosen 
such that their mass dependence is also of the form 

4 = ln^, 4 = ln^. (20.5) 

The constants fc^, kp collect the analytic contributions at order x"^ and are in- 
dependent of m and /i. By completing the logarithms at oder x"^ to a square, 
we have in effect chosen a particular running scale: the one where the coefficient 
of the term linear in InM^ vanishes. This simplifies the representation, but is 
without physical significance - the decomposition into an infrared singular part 
arising from the Goldstone bosons and a regular remainder is not unique. In 
ref. [BT|, a somewhat different representation for the analytic terms of order x^ is 



used, which involves the two parameters r^, Vp instead of k^, kp. 



21 Numerical results for quark mass dependence 

In addition to the coupling constants £1, ... ,£4 that also govern the low energy 
properties of the vrvr scattering amplitude, the ratios M^/M^ and Fj^/F contain 
the two fourth order constants /cai and kp. We expect that the contributions from 
these terms are of order M^/Mj and can just as well be dropped - unless m is 
taken much larger than in nature. We did not make an attempt at quantifying 
the uncertainties associated with these terms, because they are small compared 
to those from the coupling constants £i,...,4- The bands shown in fig. |10| 
correspond to = = 0. If we were to instead set the constants /cm, kp 
equal to zero, the boundaries would be slightly shifted, but the shifts are small 
compared to the width of the bands. 

For small values of m, the contributions of O(m^) dominate. These are de- 
termined by the two scales A3 and A4. As discussed in section |^, the infor- 
mation about the first one is meagre - the crude estimate ( [7.8|) amounts to 
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0.2 GeV < A3 < 2 GeV. For the second one, however, the value for £4 obtained 
in section |TT| yields a rather decent determination: 



A. = 1.26 ±0.14 GeV. 



(21.1) 



The two parameters A3, A4 play the same role as the coefficients Cm and Cp in the 
polynomial approximations = 2 mB [1 + Cm m) , F^^ = F {1 + Cp m) , that are 
sometimes used to perform the extrapolation of lattice data. In contrast to these 
approximations, the formulae (|20.4|) do account for the infrared singularities of 
the functions M.n-(m) and FT^{m) - to the order of the expansion in powers of m 
considered, that representation is exact. 

Consider first the ratio F^^/F, for which the poorly known scale A3 only enters 
at next-to-next-to-leading order. The upper one of the two shaded regions in 
fig. |1^ shows the behaviour of this ratio as a function of M, according to formula 
(|20.4 ). The change in occurring if M is increased from the physical value to 
Mk is of the expected size, comparable to the difference between Fk and F^. 
The curvature makes it evident that a linear extrapolation in m is meaningless. 
The essential parameter here is the scale A4 that determines the magnitude of 
the term of order . The corrections of order are small - the scale relevant 
for these is kp — 0.5 GeV. 




0.2 0.3 

M(GeV) 



Figure 10: Dependence of the ratios F^^/F and M^/M"^ on the mass of the two 
lightest quarks. The variable M is defined by M ^ = (m^ + m^) B and F is the 
value of the pion decay constant for m„ = nid = 0. The strange quark mass is 
held fixed at the physical value. 



In the case of the ratio M^/M^, on the other hand, the dominating contri- 
bution is determined by the scale A3 - the corrections of O(M^) are small also 
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in this case (the relevant scale is Aa/ — 0.6GeV). The fact that the information 
about the value of £3 is very meagre shows up through very large uncertainties. 
In particular, with A3 ~ 0.5 GeV, the ratio M^/M^ would remain very close to 
1, on the entire interval shown. Note that outside the range of values for £3 con- 
sidered in the present paper, the dependence of on the quark masses would 
necessarily exhibit strong curvature. This is illustrated with the dashed line that 
indicates the behaviour of the ratio M^/M^ for £3 = —10. According to fig. ^ 
this value corresponds to aj] — 0.24. 

The above discussion shows that brute force is not the only way to reach the 
very small values of m„ and observed in nature on the lattice. It suffices 
to equip the strange quark with the physical value of and to measure the 
dependence of the pion mass on mu,md in the region where M,r is comparable 
to Mk- Since the dependence on the quark masses is known rather accurately 
in terms of the two constants B and A3, a fit to the data based on eq. ( |20.4| ) 
should provide an extrapolation to the physical quark masses that is under good 
control. Moreover, the resulting value for A3 would be of considerable interest, 
because that scale also shows up in other contexts, in the tttc scattering lengths, for 
example. For recent lattice work in this direction, we refer to |37|. A measurement 
of the mass dependence of F^, in the same region would be useful too, because it 
would provide a check on the dispersive analysis of the scalar radius that underlies 
our determination of A4 - in view of the strong unitarity cut in the scalar form 
factor, a direct evaluation of the scalar radius on the lattice is likely more difficult. 
Chiral logarithms also occur in the quenched approximation but since 

the coefficients differ from those in the full theory, a naive comparison of the 
above formulae with quenched lattice data is not meaningful. 

In section ^ we noted that the expansion of the scattering length Oq in pow- 
ers of the quark mass contains an unusually large infrared singularity at one loop 
level. We can now complete that discussion with an evaluation of the contribu- 



tions arising at next-to-next-to-leading order, repeating the analysis of ref. |^ 
with the information about the coupling constants available now. The result is 
shown in figure |Tl|, where we indicate the behaviour of the correction factors 
-Ro, -R2, defined by 

Wn — -fLn , tin — TTTT -fto , 

° 32 7rF2 ° 16 ttF^ ' 

as a function of M^^. The reason for choosing the variable rather than M is 
that the uncertainties in A3 then affect the result less strongly. The comparison 
with figure |l^, where a larger mass range is shown, demonstrates that Rq grows 
much more rapidly with the quark masses than F^^. The effect arises from the 
chiral logarithms associated with the unitarity cut - the coefficient of the leading 
infrared singularity in aj] exceeds the one in F^^ by a factor |. Note that the 
chiral perturbation theory formulae underlying the figure are meaningful only in 
the range where the corrections are small. The shaded regions exclusively account 
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Figure 11: Size of the corrections to Weinberg's leading order predictions for 
the TTTT S-wave scattering lengths, as a function of the pion mass. 



for the uncertainties in the values of the coupling constants. In the case of Rq, 
those due to the terms of order are by no means negligible for M^. > 0.2 GeV, 
so that it matters what exactly is plotted. The curves shown in the figure are 
obtained by expressing Rq, R2 in terms of the coupling constants F, ii, . . . and 
of M^, expanding the result in powers of M^^ and truncating the series at order 
M^. On the left half of the figure, the behaviour of Rq obtained, for instance, by 
truncating the expansion in powers of M instead of the one in M^^ is practically 
the same, but on the right half, there is a substantial difference, indicating that 
the chiral perturbation series is out of control there. 

We conclude that in the case of the 1 = scattering length, a meaningful 
extrapolation of lattice data to the physical values of m„ and uid requires sig- 
nificantly smaller quark masses than in the case of or F^r. In the 1 = 2 
channel, the effects are much smaller, because this channel is exotic: The final 
state interaction is weak and repulsive. The lattice result, Oq = —0.0374 ±0.0049 
7|], corresponds to R2 = 0.82 ± 0.11. It is on the low side, but not inconsistent 



with our prediction: = -0.0444 ± 0.0010, R2 = 0.98 ± 0.02. As in the case of 
and M^^, the comparison between the lattice result and our prediction is not 
really meaningful, because that result relies on the quenched approximation. The 
evaluation of the scattering lengths to one loop in the quenched approximation 
[ |69| , [72| has shown that the infrared singularities are different from those in full 



QCD. Moreover, as pointed out by Bernard and Golterman |72|, the very method 



used to extract the infinite volume scattering lengths from finite volume observ- 
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ables is affected: In addition to the purely statistical error, the numbers in 
|n]] have a sizeable systematic error. 



22 Summary and conclusion 

1. The Roy equations determine the vrvr scattering amplitude in terms of the 
imaginary parts at intermediate energies, except for two subtraction constants: 
the S-wave scattering lengths Og, Og. At low energies, the contributions from the 
imaginary parts are small, so that the current experimental information about 
these suffices, but the one about the scattering lengths is subject to comparatively 
large uncertainties. 

2. The low energy theorems of chiral symmetry provide the missing element: 
They predict the values of the two subtraction constants. In the chiral limit, 
where the pions are massless, the S'-wave scattering lengths vanish. The breaking 
of the symmetry generated by the quark masses m„ and rud leads to nonzero val- 
ues, for as well as for Og and Og. The Gell-Mann-Oakes-Renner relation ( |1 . 1| ) 
shows that the leading term in the expansion of in powers of the quark masses 
is determined by the quark condensate and by the pion decay constant. Wein- 
berg's low energy theorems (|1.2| ) demonstrate that the same two constants also 
determine the leading term in the expansion of the scattering lengths. Ignoring 
the higher order contributions, these relations predict af] = 0.159, Og = —.0454. 

3. Chiral perturbation theory allows us to analyze the higher order terms 
of the expansion in a systematic manner. In the isospin limit, m„ = 171^ = m, 
the perturbation series of the vrvr scattering amplitude has been worked out to 
next-to-next-to-leading order (two loops). The result, in particular, specifies the 
expansion of Og and ag in powers of m up to and including 0{m^). The isospin 
breaking effects due to m„ 7^ ma have also been been studied [0. These effects 
only show up at nonleading orders of the expansion and are small - in contrast 
to the kaons or the nucleons, the pions are protected from isospin breaking. In 
the present paper, we have ignored these effects altogether. 

4. Chiral symmetry does not fully determine the higher order contributions, 
because it does not predict the values of the coupling constants occurring in 
the effective Lagrangian. There are two categories of coupling constants: terms 
that survive in the chiral limit and symmetry breaking terms proportional to a 
power of m. The former show up in the momentum dependence of the scattering 
amplitude, so that these couplings can be determined phenomenologically. For 
the couphng constants of the second category, which describe the dependence on 
the quark masses, we need to rely on sources other than vrvr scattering. 

5. The higher order terms of the expansion are dominated by those of next-to- 
leading order, which involve the coupling constants ii, £2 from the first category 
and £3, £4 from the second. We rely on the dispersive analysis of the scalar pion 
form factor to pin down the coupling constant £4. Crude theoretical estimates 
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indicate that the contributions from £3 are very small, but the uncertainties of 
these estimates dominate the error in our final result for the S'-wave scattering 
lengths. We also rely on theoretical estimates for the symmetry breaking coupling 
constants ri, . . . , ,r4 of next-to-next-to-leading order. These indicate that the 
contributions to Oq and from those constants are tiny and could just as well 
be dropped. 

6. The expansion in powers of the quark masses contains infrared singularities 
- the chiral logarithms characteristic of chiral perturbation theory. In the case 
of Qq, for instance, these singularities enhance the magnitude of the corrections 
quite substantially. The origin of the phenomenon is well understood: The final 
state interaction in the / = S'-wave generates a strong branch cut. For this 
reason, the straightforward expansion of Oq in powers of m converges only rather 
slowly. We exploit the fact that, in the subthreshold region, the expansion of the 
scattering amplitude converges much more rapidly: In our approach, the chiral 
and phenomenological representations of the scattering amplitude are matched 
there. With this method, even the one loop approximation of the chiral pertur- 
bation series yields values for the scattering lengths that are within the errors of 
our final result, which reads 

a°o = 0.220 ± 0.005 , al = -0.0444 ± 0.0010 , 

2a° - 5al = 0.663 ± 0.007 , - ^0 = 0-265 ± 0.004 . 

7. We have worked out the implications for the phase shifts of the S- and 
P-waves. As shown in figs. 0-^, chiral symmetry and the existing experimental 
information constrain these to a rather narrow range. The corresponding predic- 
tions for the scattering lengths of the D and F-waves, as well as for the effective 
ranges are listed in table 

8. Our representation of the scattering amplitude, in particular, also yields 
an accurate prediction for the phase of e'/e: The result for the phase difference 
between the two S-waves at s = M]^o reads 

S^oiMlo) - 5o(^io) = 47.7° ± 1.5° . (22.1) 

9. The mass and the width of the p-meson can be calculated to within 
remarkably small uncertainties: 

= 762.4 ± 1.8 MeV, = 145.7 ± 2.6 MeV . (22.2) 

Also, we confirm that the / = S'-wave contains a pole far away from the real 
axis, at ^/s = (470 ± 30) — i (295 ± 20) MeV. The phenomenon is related to the 
fact that chiral symmetry requires the scattering amplitude to be very small at 
threshold and then to grow with the square of the energy. 

10. The consequences for the coupling constants ii, £2, ^4, ^5 and rg of the 
effective Lagrangian were studied in detail. Our results are in good agreement 
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with previous work, but are more accurate. In particular, we have shown that 
£i and £2 are accompanied by strong infrared singularities generated by the two 
loop graphs, which shift the numerical values obtained in one loop approximation, 
quite substantially. The effective couplings relevant at one loop level are given in 
eq. ( |14.3| ). We have shown that, with these values, a very decent representation 
of the scattering amplitude is obtained by matching the Roy equations with the 
one loop approximation of chiral perturbation theory: The result can barely be 
distinguished from the representation that underlies the present paper. 

11. The results for the various quantities of interest are strongly correlated. 
We have examined the correlations in detail, not only for the threshold parameters 
and coupling constants, but also for the coefficients 61, . . . ,bQ of the polynomial 
that occurs in the chiral perturbation theory representation of the scattering 
amplitude. 

12. The resulting picture for the low energy structure of the scattering am- 
plitude is consistent with the resonance saturation hypothesis: For all of the 
effective coupling constants encountered in the two loop representation, our re- 
sults are consistent with the assumption that, once the poles and cuts due to pion 
exchange are removed, the low energy structure of the amplitude is dominated by 
the singularities due to the lightest non-Goldstone states 0. Since the splitting 
of a resonance contribution from the continuum underneath it is not unique, the 
saturation hypothesis does not lead to very sharp predictions, but it is by no 
means trivial that these are consistent with the values found, both in sign and 
magnitude. 

13. On the lattice, it is difficult to reach the small values of rriu and that 
are realized in nature. We have shown that chiral perturbation theory can be used 
to extrapolate the results obtained at comparatively large values for these masses, 
in a controlled manner. The method at the same time also allows a measurement 
of some of the coupling constants that occur in the symmetry breaking part of the 
effective Lagrangian, in particular, of £3. In quite a few cases, the uncertainties 
in our results are dominated by those from this term. 

14. We emphasize that most of our results rely on the standard picture, ac- 
cording to which the quark condensate represents the leading order parameter 
of the spontaneously broken symmetry, so that the Gell-Mann-Oakes-Renner re- 
lation holds. The crude theoretical estimates for the coupling constant £3 we 
are relying on indicate that the higher order terms in the expansion of 

are very small, so that the square of the pion mass indeed grows linearly with 
m = I {rriu + md) - a curvature only shows up if m is taken much larger than 



the physical value. In ref. p3[, £3 is instead treated as a free parameter and is 
allowed to take large values, so that the dependence of on the quark masses 
fails to be approximately linear, even in the region below the physical value of 
m. There is no prediction for the scattering lengths in that framework. 

15. Even if the quark condensate is not assumed to represent the leading order 
parameter, a strong correlation between Oq and Oq emerges, which originates in 
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the relation between these quantities and the scalar radius. The correlation is 
of interest, in particular, in connection with the analysis of K^.^ data: As shown 
in fig. 1^, the preliminary results of the E865 experiment at Brookhaven [^, 30 



yield a remarkably good determination of Oq. The outcome beautifully confirms 
the prediction (|11.2| ): The best fit to these data is obtained for = 0.218, with 
= 5.7 for 5 degrees of freedom. For a detailed discussion of the consequences 
for the value of a^, we refer to |^6|, 



16. A measurement that aims at determining the lifetime of a tt+tt^ atom 
to an accuracy of 10% is currently under way at CERN. The interference of the 
electromagnetic and strong interaction effects in the bound state and in the decay 
is now well understood, also on the basis of chiral perturbation theory The 
decay rate of the ground state can be written in the form 

r = '^a'p\a^o-al\'{l + 6), 



with p = y — M^o — \o? . The term 6 accounts for the corrections of 
order a as well as those due to m„ 7^ niii. According to ref. these effects 



increase the rate hj S = 0.058 ± 0.012, that is by about 6%. Inserting our result 
( |11.2| ) for Og — clq, we arrive at the following prediction for the lifetime : 

r = (2.90 ±0.10) ■ 10"^^ s. (22.3) 

Since the decay rate is proportional to |aQ — Oq^, the outcome of the experiment 
is expected to lead to a determination of |ag — Ogl to an accuracy of 5%, thereby 
subjecting chiral perturbation theory to a very sensitive test. 
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A Notation 

We use the following abbreviations: 



e={ -] , x={ , L = ln^, N=16n\ 

The intrinsic scales of the coupling constants of £4 are denoted by A„. In terms 
of these, the standard renormalized couplings are given by 



A^ 

lif^) = ^ In — , 71 = 3 ' 72 = 3 , 73 = -2 , 71 
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where /i is the running scale. The scales relevant in the various applications are 
different and the formulae can be simplified considerably if the coupling constants 
are normalized at the appropriate scale. For this reason, we use three different 
symbols: 

/ = In / = In — / = In — 

The first coincides with the one introduced in 0], in the framework of a one loop 
analysis, where there is no need to distinguish from £„. The quantities differ 
from the running coupling constants ^^(/u) only by a numerical factor, which it is 
convenient to remove in order to simplify the expressions. For the same reason, 
we work with the coefficients &i, . . . , fee, introduced in 0, 

h = Nh, b2 = Nb2, h = Nh, h = Nb4, b5 = N%, be = N%. 

and also stretch the coupling constants of Cq by a power of = 16 vr^, 

f„ = Ar2r;(/x), (n = l,... ,8). 

In Generalized Chiral Perturbation Theory, the coefficients bn are replaced by 
the constants a, P, Xi, . . . , A4: 

a = 1 + ^3 b, + Ab2 + Ah - Ah) - '^^n' e - ^ 

/5 = 1 + ^(62 + 4 63 -4 64) +4^' (3 h-h)-§n'e+'-f^'^ 

NX2 = 2h-l^, (A.l) 
A^'As = 65 - I &6 + i , 

/\r2x _ 4 T 14 
iV A4 — — 3 Og — 27 • 

Throughout this paper, we identify with the mass of the charged pion and 
use = 92.4 MeV 



B Polynomial part of the chiral representation 

In this appendix, we convert the representation obtained in ref. for the scat- 
tering amplitude to two loops of chiral perturbation theory into an explicit ex- 
pression for the coefficients of the polynomial C{s,t,u) defined in eq. (^I6|) . As 
a first step, we decompose that representation into three functions of a single 
variable and a polynomial, according to eq. (|2.5|) . The resulting representation 
for the functions U^{s), U^{s), U'^{s) contains linear combinations of the loop 
integrals J(s), . . . , K^i^s) introduced in ref. 0. Note that the decomposition is 
not unique: eq. ( p^.5| ) fixes the functions U\s) only up to a polynomial in s. 
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Next, we expand the loop integrals in powers of s, using the explicit expres- 
sions of ref. 1^. In terms of the dimensionless variable s = s/M^, the result 
reads: 

A-.W = -^U + ^ + lS + o(.»)^ (B.i; 

Ms) -^,{ — + 30 + 420 + ') • 

^*<'''""3072^\~3~ + 20 + 105 + '/ ' 

These representations suffice to determine the Taylor series expansions of the 
functions U\s) around s = 0, up to and including O(s^) in the case of U^{s), 
U^{s) and to 0{s^) for U^{s). 

The ambiguity mentioned above is fixed with the dispersive representation 
( |2.7| ), which requires the first few derivatives of the functions U\s) to vanish 
at s = 0. Starting with an arbitrary decomposition, for which this requirement 
need not be obeyed, we truncate the Taylor series for f/°(s), f/^(s) at order and 
the one for U^{s) at order s^. It is straightforward to check that the functions 
obtained by subtracting these terms indeed obey the dispersion relations (|2.7|) . 
Absorbing the subtractions in C{s,t,u), we may then read off the coefficients of 
this polynomial: 

Cl = 




C2 




59 62 731 63 

C3 = 7 < b3 + — + n + + 

3 nrr.A ^ 3 ^ , . ^ I ^ 105 315 



3601 64 5387 TT^ 19121\1 
^ 315 15120 ~ 7560 / f' ^ ' ' 



= J—lb 31 / 43 h 8 62 23 63 
^4 jY^4 \ 4 2520 H 420 63 63 



48 



997 64 467 _ 63829 \1 
^ 315 ^ 7560 ~ 45360/ J ' 
_ 1 f 137 bi 379 h 25 63 731 64 ^ 
^5 - jv2^6 j 1680^ 16 1680 28 ISCT ^ ^ 

269 616731 
^ 15120 ^ 18144 J ' 

_ 1 f 31 bi 4762 6563 54764 T 
- jY2^6 j~i680^ ^ 112 ~ 1680 ~ '252 420 ^ ^ 

2 



TT- 442871 
^15120 ^ 90720 J ■ 

The constants 6„ represent dimensionless combinations of coupling constants, 
introduced in ref. In the notation of appendix 0, the expressions read: 

T 7L 4£i I, ^.13 m^L' 

01 = 2c4H 

6 3 2 18 6 

4£i 56^2 ,~ 26^4 47 1 IGi^i, ^ 

-3£3^~4-5£7 + ^ + ^-^ + ^-^ + ^|, 
^ ^ 27 27 4 9 27 1296 J ' 

3 3 9 36 ^\ 9 2 3 54 J 

{ 16 , - - , ~ 166£2^94 8£~4^3177r2 17891 

— ^+^^^^+'^^-'^^-^+—-^+^-13^1' 

63^^ + ^l + ^-l-^ + eLj-^-^ + 2£4 + ^| 
236 12 9 \9 3 216j 

/ 4£i£4 , 2^2^ , 89£i , 38^2 7^4 SIItt^ 70631 ,^ 
- L ^ £2 5 5^L^ f£i , 84 , 2^4 47 ] 

H T h ^-/^ 



6 6 36 6 ^ [ 9 9 3 216 J 

+ Jf 5£i 4^2 54 177r^ 1655] 

^ 3 27 27 9 216 2592 J ' 

- 85L2 ~(7ii 107^2 6251 31£i 145^2 7 66029 
05 = h < \ } + 1 , 

72 \ 8 72 288 J 36 108 54 20736 

- ff5£i 25^2 2571 7ii 35^2 tt^ 11375 
u I- T ■{ - -I ^- > 4- f - - -I 

^24 \ 72 72 864 J ^ 108 108 27 20736' 

Note that the quark masses exclusively enter through ^ and L - the remaining 
quantities are independent thereof. The expressions involve the logarithm of M^, 
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as well as the square thereof - in the chiral limit, the coefficients 6„ diverge 
logarithmically. The coefficients of the leading infrared singularities are pure 
numbers, which are determined by the structure of the symmetry group and the 
transformation properties of the symmetry breaking part of the Hamiltonian, 
that is of the quark mass term. The scales of the logarithmic divergences, on the 
other hand, are not determined by the symmetry, but are fixed by the intrinsic 
scales Ai, . . . , A4 of the effective coupling constants of £4. 



C The corrections Ao, Ai, A2, 

The leading terms in the expansion of the scalar radius in powers of the quark 
masses are determined by £4. The next-to-leading order correction A^,, which is 
defined in ( |7. 5| ) , was calculated in on the basis of an evaluation of the scalar 
form factor to two loops. In the notation introduced above, their result reads: 



A , = - ^ L 2 + L ( - ^ - ^ £2 + 4 £4 - ^ I + r 5, + 4 + 2 



18 1 9 9 216 J 

22 ^7 ^7, 13 r, 23 869 

H £1 + 2^2 4 4 \ (C.l) 

9 24 6 72 648 ^ ' 

As discussed in section ||, the quantities Co, Ci and C2 tend to unity in the 
chiral limit. According to equation (|7.6| ), the ffist order corrections can be ex- 
pressed in terms of the scalar radius and the coupling constant £3. To work out 
the corrections of second order, it suffices to insert the relations (|B.2|) , ( p.3|) in 
the definition ( |6.2[ ) of Ci and C2 and read off the coefficients of the terms of order 

The result reads 

71L2 40£i 804 54 ~ 53931 

Ai = + L{ ^+4^4 + ) 

12 \ 9 9 2 ^ 315 J 

+ fa + 4f3 - 4^4 - 2f52 - 44 + 42 

1826 4 3118 4 794 144 4 521 . 24221 ,^ , 

H \ \ TT^H , C.2) 

315 315 21 35 252 3024 ' ^ ' 

175 ~ \ ~ 148 4 ~n 93111 

A2 = + L{-lQli + 4 + 64 + ) 

18 \ 9 630 J 

i"^ ~ ~ 

- 4 + 44 - 44 - 2 + y + 44 + 4' 

556 4 43724 1044 1254 2939 TT^ 43109 

+ + + + . 

63 315 35 21 1260 5040 

According to ( |6.6| ), the analogous correction in the low energy theorem for Cq 

can be expressed in terms of these: 

Ao = ^(12Ai-5A2). (C.3) 
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D Phenomenological representation 

In the present appendix, we convert the low energy representation of the scat- 
tering amphtude constructed in ref. ||^ into the form given in section |^. That 
representation consists of a sum of two terms: 

A{s, t, u) = A{s, t, u)sp + A{s, t, u)d , (D.l) 

The first describes the contributions generated by the imaginary parts of the S- 
and P-waves below y/s2 = 2 GeV, while the background amplitude A{s, t, u)d 
collects those from the higher partial waves and higher energies. 

The explicit expression for the first term involves three functions of a single 
variable: 

A{s, t, u),p = 32n {lW\s) + |(s - u)W\t) + |(s - t)W\u) 

+ ^W'{t) + ^W\u) - lW\s)} . (D.2) 

The functions W'^{s), W^{s), W'^{s) are determined by the imaginary parts of 
the S- and P-waves and by the two subtraction constants Oq, a^: 



W'^is) = + ^ ^ / ds' 



AMI JAMi s'(s'-4M2)(s'-s) ' 

W\s)='^r ds' "^f^, (D.3) 



4M2 TT JiMi s'{s' -AM^){s' - s) 

These functions are closely related to the quantities W\s) introduced in section 
1^, but there are two differences: The subtractions are not the same and the range 
of integration differs. 

To compare W^{s) with W^{s), we consider the function 



w [s 



71 JiMl S"^{s' — s) 



which is intermediate between the two: It differs from W^{s) only in the range 
of integration and from W^{s) only by a subtraction polynomial. The latter may 
be expressed in terms of the following moments of the imaginary part: 

n 1 /■■^2 dslmtnis) 
Jn = - —nr, VTT^' n = 0,l,2. p.4 

The explicit relation between W^{s) and w^{s) reads 

W\s) = w%s) + ^ + sis- AMI) J'o +s'is- AMI) ^? + 4M^V° . (D.5) 

TT 
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The difference W^{s) — on the other hand, is given by an integral over 

the region s > S2, which merely generates contributions of 0{p^), so that 

w\s) = W'^is) + 0(/) . 

To the accuracy to which the two loop representation holds, we may thus replace 
the term w^{s) on the r.h.s of eq. ( |D.5| ) by W^{s). 

This shows that, up to terms of 0{p^), the functions W^{s) and W^{s) only 
differ by a polynomial whose coefficients are determined by the moments J°. 
The same reasoning also applies to the components with / = 1,2. The relevant 
moments are given by 

,^1 dslvcit\{s) r-^ dsluitl{s) 

" vr 74A/2 _ 4M2) ' " vr J^Mi - 4M2) ■ ^ ■ ^ 

The net result amounts to a representation for A{s,t,u)gp in terms of the func- 
tions W^{s) and a set of polnomials involving the above moments: 

W\s) = W\s)+ ^,+s{s- AMI) 4 + s\s - AMI) Jl + 4M^^V° + 0(/) , 
W\s) = W\s) -sJl- s^Jl + 0(/) , (D.7) 
W\s) = W\s)+ ^ + sis- AMI) Jo + s'is - AMI) Ji + AM^jj + 0(/) . 



E Moments of the background amplitude 

We now turn to the second part in the decomposition ( p.l| ). The chiral represen- 
tation shows that the infrared singularities contained in A{s, t, u)d start manifest- 
ing themselves only at higher orders. Up to and including O(p^), the background 
amplitude is a crossing symmetric polynomial of the momenta: 

A{s,t,u)d = P{s,t,u) + 0{p'') (E.l) 

P(S, t,u) = P1+P2S+ P3 +P4{t- uf + J95 + P6 s{t - u)"^ . 

The coefficients pi, . . . , pq may be calculated by expanding the dispersion inte- 
grals in powers of the momenta. For a detailed discussion, we refer to appendix B 
of ref. [^]. By construction, A{s,t,u)d does not contribute to the S'-wave scatter- 
ing lengths. This condition fixes pi and p2 in terms of the remaining coefficients: 

Pi = -16M>4, P2 = AM^{-p3+Pi-AM^p,). (E.2) 

The explicit expressions for these read 

Ps = ^ « - 9/0^ - l!) + ^ (-8 /? -21/^ + 11 + 12 H) 
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J94 = 8 vr {I^ + J2) + 16 vr (/^ + if) (E.3) 
P5 = ^(8/° + 9/|-ll/2-6if) 
P6 = 4 7r(/^ -3/i2 + 2if). 

The moments 1^ represent integrals over the imaginary parts at t = (total 
cross sections), except that the contributions from the S- and P-waves below 
^/s2 = 2 GeV are removed. In terms of the imaginary parts of the partial waves, 
the explicit expression reads 0: 

ri v (2/ + l) p, Imti(.) 
^ g (2/ + 1) foo^^ Iratiis) 



£=0 



TV Js2 S"+2(S - 4M2) 



The term H represents an analogous integral over the derivatives of the imaginary 
parts with respect to t at t = 0. Since the S-wave contributions are independent 
oft, they drop out. Moreover, on account of crossing symmetry, the contributions 
with 1 = 1 may be expressed in terms of those with / = 0, 2: 

,r, Tc JiMi 9 s'^ s — 4Mi 



1=2 



There is an important difference between the moments relevant for the back- 
ground amplitude and those associated with the 5"- and P-waves: While and 
H remain finite when the quark masses are sent to zero, 

/o^ = 0(l), A' = 0(1), /| = 0(1), H = 0{1), 

the S- and P-wave moments with n > 1 explode in that limit: 

4 = 0(1) , J[ = 0{M-') , Ji = 0{M~') . 

The phenomenon arises from the manner in which we have chosen to decompose 
the amplitude into a contribution from the S- and P-waves, described by the 
functions W\s), and a polynomial. These functions develop an infrared singu- 
larity in the chiral limit, which cancels the one occurring in the polynomial - 
the full scattering amplitude approaches a finite limit when the quark masses 
are turned off. In fact, the problem does not occur in the original form of the 
decomposition, based on the functions W^{s): these do have a decent chiral limit. 

The same singularities also show up in the chiral representation of the am- 
plitude: As we have normalized the unitarity correction by subtracting the dis- 
persion integrals at s = 0, they contain a quadratic infrared divergence in the 
chiral limit. Indeed, the relations (p.2| ) show that the coefficients C5 and Cg con- 
tain contributions that are inversely proportional to and precisely cancel this 
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divergence. The above choice of the decomposition has the advantage that it is 
scale independent and leads to a simple form of the matching conditions. The 
backside of the coin is that the two pieces do not have a smooth chiral limit. We 
could modify the normalization of the unitarity correction in such a way that it 
remains finite in the chiral limit, at the price of introducing an arbitrary scale to 
normalize the logarithmic infrared singularities. There is no gain in doing that, 
however: Anyway, only the sum of the polynomial and the unitarity correction 
is relevant, so that different decompositions lead to identical results. We stick to 
the one above. 

Finally, we add up the two parts of the amplitude. The result takes the form 
of eq. ( p.2|) , where P(s, t, u) represents the sum of the two polynomials associated 
with the two parts. The coefficients pi, . . . , involve a linear combination of the 
various moments introduced above. In fact, up to terms of O(p^), the result may 
be expressed in terms of the combinations 

^ n ~ ~^ '^n ^ ^ n ~ ~^ , I ^ = -\- . (E.6) 

The contributions from the S- and P-wave moments precisely represent the 
pieces needed to complete the sum over the angular momenta - the factor of 3 in 
front of the P-wave moments accounts for the weight 2£ + 1 that occurs in the 
definition ( [E.4|) of the moments J^. The result is independent of the energy S2 
used in the decomposition ( |D.1| ) of the amplitude and is given in eqs. ( |3l^ ) , ( |3.5| ) . 

The moments are readily evaluated with the information given in ref. 0. 
Since the angular momentum barrier suppresses the higher partial waves near 
threshold, the terms Jj, I2 and /| are negligibly small. The contributions from 
the S- and P-waves depend on the values of the two S'-wave scattering lengths 
and Oq. In the narrow range of interest here, this dependence is well described 
by a quadratic interpolation of the form 

I = Uo + Ui Aa° + U2 Aoq + M3 (Aoq)^ + M4 Aog Aog + M5 (Aag)^ , 

with Aoq = — 0.225, Aoq = Og + 0.03706. The numerical values of the coeffi- 
cients are listed in table ^. 





1 = 


/ = 1 


1 = 2 




^ 


n 

GeV"^ 


-' 2 
GeV"^ 


n 

GeV"^ 


GeV"^ 


GeV"^ 


^ 

GeV"^ 


n 

GeV"^ 


-' 2 
GeV"^ 


Ul 
U2 
U3 
U4 

U5 


9.44 

58.5 
-75.2 

82.4 
-232 

322 


66.7 
507 
-462 
902 
-1920 
2280 


609 
4980 
-3300 
9800 
-16200 
16900 


1.90 
1.97 
-7.99 
-1.09 
-2.66 
17.4 


3.92 

6.35 
-25.2 

-2.1 
-18.3 

72.6 


10.6 
25.6 
-99.3 
-2.84 
-110 
349 


.469 
1.56 
-15.1 
-.295 
-14.5 
125 


2.61 
5.84 
-97.8 
-2.65 
-54.3 
932 


21.4 
33.9 
-884 
-23.6 
-317 
9190 



Table 3: Moments of the background amplitude. 
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The quantity H exclusively receives contributions from the partial waves with 
i > 2. As discussed in detail in 0, H is dominated by the contribution from 
the lowest spin 2 resonance, which is independent of Og, a^. The behaviour of 
the integrand in the threshold region does depend on the two S'-wave scattering 
lengths, because these determine the threshold parameters of the higher partial 
waves, but since the contributions from that region are very small, we ignore the 
dependence on aj], Oq and use the value given in 

i^ = 0.32GeV"^ (E.7) 

F Error analysis and correlations 

The matching conditions and the two loop representation for the scalar radius 
determine the values of the constants 

as functions of the parameters occurring in our input, 

y = {{r'^)s, ^3,^i,^2,r3,r4,r52,6'o,6'i} . 

The last two of these characterize the experimental input used when solving the 
Roy equations. Strictly speaking, that part of our input involves three functions, 
the imaginary parts of the S- and P-waves above 0.8 GeV, but in practice, the 
solution of the matching conditions is sensitive only to two parameters: the values 
of the phases Sq{s) and Sl{s) at y/s = 0.8 GeV, which we denote by 9q and 6*1, 
respectively. In ref. the uncertainties in these parameters are estimated at 
^0 = 82.3° ± 3.4°, 01 = 108.9° ± 2°. 

The uncertainties in the input give rise both to uncertainties in the individual 
components of x and to correlations among these. We describe the correlations 
in terms of a Gaussian distribution: The probability for x to be contained in the 
volume element dx is represented as 

dP = N exp{-lQ) dx, (F.l) 
Q = J2Cab Ax" Ax" , Ax'' = x''- {x'') . 

ab 

The uncertainties in our results are characterized by the coefficients Cab of the 
quadratic form in the exponential. In the Gaussian approximation, the matrix C 
is given by the inverse of the correlation matrix K, 

J^ab ^ (^^^a^^b'^ ^ C-K=l, (F.2) 
so that the error analysis boils down to an evaluation of the matrix K. 
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In the small region of interest, the response to a change of the input variables is 
approximately linear. We denote the central values of the input parameters by Uc 
and set = (1 + //*), with (jf) = 0. Linearity then implies that the mean value 
of x'^ coincides with the solution of the matching conditions that corresponds 
to our central set of input parameters, and that the correlation matrix K can 
be expressed in terms of the matrix {rj'^r]^). We treat the input variables as 
statistically independent, so that this matrix is diagonal, 

{v'v') = S'\(^^Y■ (F.3) 

The result of the dispersive evaluation, (r^)^ = 0.61±0.04fm^, implies that the 
mean square deviation in the variable 1]^ is given by cxi = 0.04/0.61. We interpret 
the estimate £3 = 2.9 ± 2.4 in the same manner: (T2 = 2.4/2.9. Concerning the 
variables ri, . . . ,r4 and rs2, we assume that all values in the interval from to 
twice the value obtained from resonance saturation are equally likely, so that, 
for n = 3, . . . , 7, the mean square deviation becomes cr„ = 1/V3. Finally, the 
uncertainties in the input parameters 6q and 61 amount to erg = 3.4/82.3 and 
CTg = 2/108.9, respectively. 

We also need an estimate for the sensitivity of our results to the value of the 
scale II used when applying the resonance estimates. For the mean values, we use 
/i = Mp. To estimate the uncertainties due to that choice, we evaluate the shift 
occurring in the quantity of interest if the coupling constants r„ are held fixed, but 
the scale is replaced by /i = 1 GeV, repeat the calculation for = 0.5 GeV, and 
take the mean square of the two shifts. Likewise, the correponding contribution to 
the correlation matrix is the average of the two matrices associated with those two 
shifts. Alternatively we could assume that all values of /i in the interval between 
0.5 and 1 GeV are equally likely. The uncertainties then become somewhat 
smaller, but it suffices to slightly stretch the interval for the outcome to be nearly 
the same. 





Aa% 


Aal 


All 


A £2 


A £4 


Afs 


Afe 




2.0-10"^ 


3.2-10-^ 


1.9-10-^ 


-1.7-10-^ 


4.2-10-^ 


3.2-10-^ 


3.3-10-5 


Aal 




9.7-10-^ 


1.6-10-^ 


-1.2-10-^ 


-4.2- 10"^ 


-2.2-10-^ 


-2.2-10-5 


Ail 






3.5-10-1 


-3.3-10-2 


6.7-10-2 


-5.4-10-1 


-3.7-10-2 


AI2 








1.2-10-2 


-7.2-10-3 


1.1-10-2 


-4.6-10-3 


AU 










4.8-10-2 


-9.1-10-2 


-2.2-10-3 


Afs 












1.1 


9.2-10-2 


Afe 














1.1-10-2 



Table 4: Numerical elements of the correlation matrix (F.2). 



The elements of the resulting correlation matrix are listed in table ^ The off- 
diagonal elements are of interest only if their numerical values are comparable 
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to the product of the square roots of the corresponding diagonal entries - the 
numbers hsted are significant only insofar as this condition is met (those below 
the diagonal are omitted - the matrix is symmetric). The errors quoted in the 
various rows of table are the square roots of the corresponding contributions to 
the diagonal elements of the correlation matrix. 

If all variables except Oq and Og are integrated out, the distribution reduces 
to a Gaussian in these two variables: 

dp = n exp(-i q) daQ dal, q = Cu (Aag)^ + 2 Cu Aa^Aal + C22 (Aoq)^ , 

where the 2x2 matrix c is the inverse of the submatrix of K that collects the cor- 
relations of Oq and Oq. The result is illustrated in fig. 0: The small ellipse shows 
the standard 68% confidence limit, that is the contour where q = 1- For a careful 
analysis of the errors and correlations associated with the various numerical eval- 
uations found in the literature, we refer to ^M, where the corresponding error 



ellipses are also shown. Note that the radiative corrections in the value of F^^ are 
often not accounted for. At the accuracy under discussion, these matter, as they 
increase the results for the 5- wave scattering lengths by about two percent. 
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